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ABSTRACT 

Ultraviolet non-ionizing continuum and mid-IR emission constitute the basis of two widely used star for- 
mation indicators at intermediate and high redshifts. We study 2430 galaxies with z< 1.4 in the Extended 
Groth Strip with deep MIPS 24 /im observations from FIDEL, spectroscopy from DEEP2, and UV, optical, 
and near-IR photometry from AEGIS. The data are coupled with dust-reddened stellar population models and 
Bayesian SED fitting to estimate dust-corrected SFRs. In order to probe the dust heating from stellar popu- 
lations of various ages, the derived SFRs were averaged over various timescales-from 100 Myr for "current" 
SFR (corresponding to young stars) to 1-3 Gyr for long-timescale SFRs (corresponding to the light-weighted 
age of the dominant stellar populations). These SED-based UV/optical SFRs are compared to total infrared 
luminosities extrapolated from 24 /im observations, corresponding to 10-18 /im rest frame. The total IR lu- 
minosities are in the range of normal star forming galaxies and LIRGs (IO'^-IO'^Lq). We show that the IR 
luminosity can be estimated from the UV and optical photometry to within a factor of two, implying that most 
z < 1 .4 galaxies are not optically thick. We find that for the blue, actively star forming galaxies the correla- 
tion between the IR luminosity and the UV/optical SFR shows a decrease in scatter when going from shorter 
to longer SFR-averaging timescales. We interpret this as the greater role of intermediate age stellar popula- 
tions in heating the dust than what is typically assumed. Equivalently, we observe that the IR luminosity is 
better correlated with dust-corrected optical luminosity than with dust-corrected UV light. We find that this 
holds over the entire redshift range. Many so-called green valley galaxies are simply dust-obscured actively 
star-forming galaxies. However, there exist 24 /im-detected galaxies, some with Lir > 1O"L0, yet with little 
current star formation. For them a reasonable amount of dust absorption of stellar light (but presumably higher 
than in nearby early-type galaxies) is sufficient to produce the observed levels of IR, which includes a large 
contribution from intermediate and old stellar populations. In our sample, which contains very few ULIRGs, 
optical and X-ray AGNs do not contribute on average more than ~ 50% to the mid-IR luminosity, and we see 
no evidence for a large population of "IR excess" galaxies. 

Subject headings: galaxies: evolution — galaxies: fundamental parameters — infrared: galaxies — ^ultraviolet: 
galaxies — surveys — galaxies: active 
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The total infrared (IR) luminosity, either al one or in combi- 
nation with the ultraviolet (UV) luminosity (iHeckman et alj 
fT998h . is increasingly being considered a reUable star for- 
matio n (SF) indicator f or normal, dusty star-forming galax- 
ies (Ke wlev et al.ll2002h . This is especially the case since the 
more traditional SF'^ indicators, such as the UV continuum 
and nebular line flux, re quire somewhat substantial correc- 
tions for dust extinction (iKennicutlll 19981) . The m/d-infrared 
luminosity has recently been suggested as a tracer of star 
formation (Roussel et al. 2001; Forster Schrei ber et al. 2004i 
Wu et al. 2005; Alonso-Herr ero et all 120061: ICalzetti et al J 
2007;,Rieke et al..,200 9), potentially serving as an alternative 
to the far IR, which is more difficult to obtain. The mid IR 
has received particular attention in intermediate and high red- 
shift studies, largely driven by the sensitivity of Spitzer MIPS 
observations, which with its 24 /im detector readily observes 
normal star form ing galaxies (Lir ~ lO'^L©) out to z ~ 1 (e.g., 
iLe Floc'h et al.l 12005) and luminous and ult ra-luminous IR 
galaxies (LI RGs, ULIRG s) out to z ^ 2 (e.g.. lPapovich et al] 
I2OO6; Reddy ltai]|2006l) . 

The validity of using the IR as a SF indicator at interme- 
diate redshifts depends critically on the assumption that the 
IR flux is tightly correlated with young stellar populations for 

'^ SF will be used to designate "star formation" or "star forming", depend- 
ing on the context. 



2 



SALIM ET AL. 



typical field galaxies in deep surveys. While one expects dust- 
reprocessed emission from both young and old stars to con- 
tribute to the IR, the question of a dominant source is less 
straightforward. The source of the far-IR emission in nearby 
star forming galaxies has been a subject of debate predating 
the launch of Spitzer Space Telescope. That the majority of IR 
heating is due to young populations, i.e., hot stars located in 
compact star forming regions, has been initially suggested by 
the similarity be tween Ha and far-IR structures within nearby 
galaxies (e.g., Dever eux et al.|[T997{ !). Studies utilizing better 
resolution from Spitzer to some degree confirmed t hese earlier 
findings and extended them down t o 70 and 24 fim dHinz et al.l 
l2004t iPerez-Gonzalez et al.ll2006l) . On the other hand, the 
claims for a more significant role of older stellar populations 
in the far IR, which heat the dust through a diffuse interstel- 
lar radiation field, were initiall y based on the modeling of 
IWalterbos & Greenawaltl (Il996h . who successfully predicted 
IRAS 60 and 100 fxm fluxes using dust models and assuming 
that Z?-band light (from intermediate age stars; ^ 1 Gyr) traces 
the general interste llar rad iation field. While it is now gener- 
ally accepted (e.g.. Ida Cun ha et al. 2008) that the interstellar 
radia tion field can be a significant heating source for the far 
IR, iBoselli et al.l (1200 Ih suggested that this may be true for 
the mid IR as well. They found that 6.75 and 15 /im emis- 
sion measured by ISO correlates better with far-IR luminosity 
than with either Ha or UV dust-corrected luminosity. More 
recently, the case for the interstellar radia tion field producing 
the 8 /im PAH emission has been made bv lBendo et al.l (l2008b 
who find a good correlation with 16 iim emission. On the 
other hand. Dfa z-Santos et al.l (12008 ^) find that 8 /im emission 
from HII regions in local LIRGs follows Paa emission from 
young stars when metallicity and age are fixed. However, un- 
like the emission at 8 /im, the general consensus for mid-IR 
continuum at 24 /im is that it is dominated by emission fro m 
star-forming regions (ICalzetti et al]l2007tlRiekeetani2009h . 

The goal of this study is to explore the use of mid-IR lu- 
minosity (specifically in 10-18 /imrest-frame range) as a SF 
indicator This wavelength range falls inbetween the 8 /im 
IRAC and 24 /im MIPS bands, where there are no direct con- 
strains from Spitzer studies of nearby galaxies. Also, our 
sample of 24 /im-detected galaxies at 0.2 < z < 1 .4 is gen- 
erally more luminous than the samples studied locally (such 
as SINGS). We base our approach on the comparison of the 
level of correlation between total IR luminosities (extrapo- 
lated from MIPS 24 /im observations) and UV/optical dust- 
corrected star formation rates (SFRs). These SFRs come from 
UV/optical SED fitting, which allows us to construct SFRs av- 
eraged over various timescales, from 0. 1 to several Gyr. SFRs 
averaged over various timescales correspond to dust-corrected 
luminosities coming from stellar populations ranging in age. 
We perform the comparison for various subsamples, specifi- 
cally for blue actively star-forming galaxies and red quiescent 
ones. Finding the age of the stellar population that best corre- 
lates with IR luminosity could indicate the stellar population 
responsible for dust heating at 10-18 /im. In ^we present the 
multiwavelength data sets used in this study. In ^we derive 
SFRs from UV/optical SED fitting, and in ^we derive IR lu- 
minosities from 24 /im observations. The results of the com- 
parison of UV/optical SFRs and IR luminosities of blue star 
forming galaxies are presented in ^ while red (dusty or qui- 
escent) galaxies and AGN candidates are analyzed in ^ In 
this paper we use a f7„, = 0.3, f^A = 0.7, Hq = 70kms~' Mpc~' 
cosmology. 



2. DATA 

In this study we use various data sets matched to the DEEP2 
redshift survey. Redshifts and UV, optical, and /T-band pho- 
tometry are part of the All-Wav elength Extended G roth Strip 
International Survey (AEGIS, iDavis et alj |2007|) . AEGIS 
combines observations from a number of ground-based and 
space observatories.'** The DEEP2 sample is /?-band selected, 
and we maintain this selection by keeping all objects even if 
they are not matched with certain bands. In most of the paper 
we study the subset of this optical sample that is detected at 
24 /im. Therefore, one has a combination of /?-band and 24 
/im selections. 24 /im data come primarily from the Far In- 
frared Deep Legacy (FIDEL) survey. Main properties of the 
data sets are given in Table [T] 

2.1. DEEP2 redshifts 

The core data set to which we match all other data is the 
DEEP2 redshift survey of the Extended Groth S trip (EGS) 
one of the four fields of the full DEEP2 survey (iDavis et alj 
2003; Faber et al. 2009). DEEP2 EGS spectra form flie basis 
of the AEGIS survey. They were obtained with the DEIMOS 
spectrograph on Keck II and cover a wavelength range 6400- 
9100 A with 1.4 A resolution. We use the 2007 version 
of the redshift catalog containing 16087 redshifts, of which 
10743 are considered secure (quality flag 3 or 4 ), represent- 
ing a 13% increase over the catalog described in lDavis et alj 
(l2007h . Galaxies were optically selected to be brighter than 
7J = 24.1''^, with a known selection function, resulting in a 
redshift distribution with a mean redshift of 0.7 and extending 
up to z ~ 1 .4 (£aber et al. 2009). 

2.2. GALEX UV photometry 

GALEX dMartin et al.ll2005l; iMorrissev et"ani2007l) imaged 
the central portion of the EGS with a single 1.2° diame- 
ter pointing (Figure [U. The exposure time was 237 ks in 
the neai--UV (NUV) and 120 ks in the fai--UV (FUV) band, 
which makes it the deepest GALEX single pointing to date. 
Data are taken from public release GR3. While GALEX ob- 
serves in FUV and NUV simultaneously, 1/2 of the FUV ex- 
pos ure time was lost due to anomaly with the FUV detec- 
tor (iMorrissey et al.ll2007h . The GALEX pipelin e produces 
catalogs using SExtractor (Be rtin & Arnouts 1996) aperture 
photometry. While adequate for more shallow, resolved 
images, such photometry suffers from severe blending and 
source confusion in the deep EGS images (GALE X resolution 
is 4 — 5", while astrometry is good to 0."5, Mor rissev et alj 
|2007|) . To remedy this problem we perform PSF sour ce ex- 
tracti on, which is less sensitive to blending (Zamoiski ejjlj 
lIOOTh . We use the custom-built PSF extraction software EM 
Photometry, developed by D. Aymeric, A. Llebaria and S. 
Arnouts, which uses the ex pectation-maximisation (EM) al- 
gorithm of iGuillaume et al] ((2006). EM Photometry extracts 
GALEX fluxes based on optical prior coordinates. While suc- 
cessfully dealing with blending, the resulting fluxes for a 
given set of objects will to some extent depend on the depth 
(i.e., the number of objects) in the prior catalog, with pho- 
tometric bias being especially pronounced for intrinsically 
fainter objects. In particular, having too many faint optical 
priors (fainter than the equivalent GALEX limit in the absence 

Please refer to |http : //aegis ■ ucolick . org| for more informa- 
tion on AEGIS, including the footprint of various data sets. 
" Magnitudes are given in AB system throughout. 
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Fig. 1 . — Areal map of the sample. Our optical sample consists of DEEP2 
galaxies (gray region) having secure spectroscopic redshifts and lying in the 
intersection of GALEX (circle) and CFHTLS (rectangle) areas (dark gray, 
0.31 sq. deg). Other surveys used in the SED fitting (i(-band MMT and K- 
band Palomar) cover the dark gray area almost entirely and their footprints 
are not shown here. Note that a galaxy is kept in the optical sample even 
if it is not detected in all bands, as long as it hes in the sample area (dark 
gray region). Thus, the optical sample is only S-band selected. We derive 
SED fitting parameters for the entire optical sample, but then study a subset 
of objects that are detected at 24 /^m. The 24 /^m data (footprint not shown) 
cover the dark gray region fully. 

of blending) will result in the splitting of the UV flux of one 
object among multiple sources, most of which are not actually 
detectable in GALEX images. We minimize this bias by us- 
ing the list of optical priors based on M-band (band closest to 
NUV) photometry from the CFHT Legacy Survey 02.3l l and 
choosing a u limit of 25.5, at which the majority of optical 
objects still have real counterparts in the NUV image. After 
extracting the NUV fluxes using the M-band prior catalog, we 
find that genuine detections mostly have NUV< 26.5, which 
we adopt as a cut for the NUV catalog, and which roughly 
corresponds to 3cr limit. As a cross-check we also per- 
form source d etection and PSF extraction using DAOPHOT 
(IStetsonlll987h . i.e., without positional constraints on source 
detections. Comparing the results from DAOPHOT and EM 
Photometry for relatively isolated sources, we find a good 
agreement for NUV< 24, and a gradually increasing differ- 
ence at fainter magnitudes, up to 0.21 mag at the catalog limit 
(DAOPHOT photometry being brighter). The difference can 
likely be attributed to unresolved detections in DAOPHOT, 
and it also represents the upper limit on the above-discussed 
bias introduced by forcing prior extractions. With an NUV 
catalog in hand, we repeat the procedure to obtain FUV pho- 
tometry, now using NUV= 26.5 to set the cut on the prior cat- 
alog and adopting a FUV= 26.5 cut for the final FUV cata- 
log (roughly a 3fT limit). We estimate the bias at the faint 
end to be below 0.13 mag. The FUV and NUV catalogs are 
matched to the CFHTLS catalog b y con struction, which is in 
turn matched to DEEP2 positions ( ^2.3l l. Of CFHTLS sources 
matched to the full DEEP2 redshift catalog, 22% have fluxes 
in FUV caalog and 59% in NUV catalog. RMS calibration 
errors of 0.0 52 and 0.026 mag are a dopted for FUV and NUV 
respectively dMorrissev et al.ll2007h . 



2.3. CFHT Legacy Survey optical photometry 

The EGS represents one of the four deep fields targeted 
by the CFHT Legacy Survey (CFHTLS). The central region 
of the EGS is observed with the MegaPrime/MegaCam im- 
ager/detector in a single pointing covering a 1° x 1° field 
of view (Figure [U in five optical bands {u*g'r'i'z')- The 
limiting magnitudes corresponding to 80% completeness are 
27.2, 27.5, 27.2, 27.0, 26.0, respectively, We use band- 
merged catalogs (publicly available version 2008A) based on 
/-band detections, with aperture photometry measured from 
/-band derived apertures. Matching to the DEEP2 redshift 
catalog is performed using a 0."3 search radius. Astromet- 
ric zero points coincide to within 0."02, and the 1-D coor- 
dinate scatter between the two catalogs is 0."08, i.e., both 
catalogs have very accurate astrometry. There are no mul- 
tiple matching candidates. Of 9923 DEEP2 objects (from 
the full redshift catalog not restricted to good quality red- 
shifts) that He within CFHTLS coverage, 9056 (91%) are 
matched. Based on the scatter of the comparison of the 
bright end with SDSS, we adopt RMS calibration errors of 
(0.04,0.025,0.025,0.025,0.025,0.035) mag for {u*g'r'i'z'). 

2.4. MMT u-band photometry 

In addition to M*-band data from CFHTL S, we also use 
M'-band photometry obtained with MegaCam (iMcLeod et alJ 
I2OO6I) on the MMT. These data extend across nearly the en- 
tire length of the EGS, with 24 overlapping fields each cov- 
ering 0.4° X 0.4°. The Sa limiting magnitude varies between 
26.3 and 27.0. Matching to the DEEP2 redshift catalog was 
performed using a 0"4 search radius after applying a 0."18 
offset in declination to bring the coordinate system of MMT 
data (based on USNO-Bl) into agreement with the SDSS sys- 
tem used in DEEP2. The 1-D coordinate scatter between the 
two catalogs is 0"12. Of 15283 DEEP2 objects (from the fuU 
redshift catalog) within MMT coverage, 10965 have a match 
(72%), with a handful of multiple match candidates, in which 
case the object with brighter u' is selected. 

Since we have M-band photometry from CFHTLS as well, 
we can compared them. The scatter between the MMT 
and CFHTLS M-band magnitudes does not increase with 
DEEP2 matching separation, indicating that the matches are 
real throughout the search radius. However, there is a 0.08 
mag overall offset between two magnitudes in the sense that 
CFHTLS M is fainter than MMT m. At the bright end we can 
compare these magnitudes to SDSS. MMT m matches SDSS 
very well, while CFHTLS u is again fainter, but by 0.05 mag 
(both MMT and CFHTLS photometry was first transformed 
to SDSS system). The offsets between MMT and CFHTLS m 
do not show an obvious color dependence. We correct these 
offsets in the SED fitting. We adopt a calibration RMS error 
of 0.04 mag for MMT photometry. 

2.5. Palomar K-band photometry 

The reddest photometry band that we use in the SED 
fitting comes from the P alomar /T-band survey of DEEP2 
dBundv et al.ll2006l l2008h . Including redder bands (such as 
IRAC 3.6 and 4.5 jjm) would not place additional constraints 
on SFRs, which are the main focus of this paper The EGS is 
almost fully covered with thirty-five 8. '6 x 8f6 WIRC frames, 
down to a 21.7-22.5 mag limit at 80% completeness. We use 

http://wwwl.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/community/CFHTLS- 
SG/docs/cfhtls.html 
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MAG_AUTO fluxes from thelBundy etaTI ( |2006|) SExtractor 
catalog, and their 1"! matching to bEEP2. Of 16087 DEEP2 
objects from the full redshift catalog, most of which are within 
K survey coverage, 10398 (65%) have a match. 

2.6. MIPS Spitzer 24 /im photometry 

In addition to UV through near-IR data that are used for 
the SED fitting, we use 24 /im observations to estimate IR 
luminosities. The 24 /im data were obtained with MIPS on 
Spitzer as part of the FIDEL survey. FIDEL observed EGS 
and ECDF-S fields with MIPS at 24 and 70 ^im to depths of 
30 /iJy and 3 mJy, respectively. These depths approach those 
of GOODS yet cover a larger area. In EGS, these data are five 
times deeper than the previous data described in Davi s et all 
(|2007) (which are co-added to FIDEL data). We extract PSF 
fluxes from 24 ^m images using DA OPHOT (MIPS has 6" 
resolution at 24 y^m: lRieke et al.l2004l) . We then match 24 /im 
sources having S/N> 3 (correspo nding to 10-1 6 /iJy) to the 
CFH12K photometry catalog (Coi l et al.1 [20041) using a l."5 
matching circle. This search radius is appropriate for bright 
24 fim sources, which have a 1-D astrometry precision of 0."5. 
However, fainter sources have poorer astrometry, so we sub- 
ject sources that initially had no match (39% of total) to a 
larger 3" radius search, recovering some 60% of them. In 
cases of multiple optical candidates (4% of cases), we pick 
the one that has the /-band to 24 /im flux ratio that is at least 
two times more likely than that of other candidates, where the 
probability is based on the flux ratio distribution of unique 
matches. This allows us to resolve 30% of multiple matches. 
We consider the remaining multiple matches to be a blend 
of more than one optical source and exclude them from the 
catalog of matched sources and from further analysis. Al- 
together, an optical match is determined for 74% of 24 /im 
sources within the optical coverage. The unmatched 24 /im 
sources are either blends or are presumably fainter than the 
R = 24.5 limit of CFH12K photometry catalog. Similar de- 
tection rates (for sim i lar R l imits) w ere found in CDF-S b y 
iPerez-Gonzalez et alJ (l2005h and by iLe Floc'h et aP (l2005h . 
70% and 60%, respectively. In the opposite direction, of 
DEEP2 objects from the fufl redshift catalog, 6581 (41%) are 
detected at 24 /im. We decide to match 24 /i m data directly 
to optical instead of using IRAC photometry dBarmbv et alj 
12008.) as an intermediate step, because IRAC coverage of EGS 
is not as extensive. As a test, for areas with IRAC coverage 
we run matching via IRAC and find that in 98% of cases we 
obtain the same optical match as with direct 24 /im to optical 
matching. 

2.7. Other data and data products 

In addition to redshifts, DEEP2 spectra provide emission 
line fluxes which we use to select narrow-line A GNs. Deriva- 
tion of fluxes is described in lWeiner et al.l ( l2007h . We also use 
Chandra X-ray detections from AEGIS-X DR2 to select X- 
ray AGN. Details of the X-ray data, cat alog construction a nd 
matching to optical sources are given in lLaird et al.l ( 120091) . 

3. UV/OPTICAL SFRS FROM SED FITTING 

The sample used in SED fitting consists of DEEP2 galax- 
ies with secure redshifts and spectra classified as galaxies. A 
small number of galaxies fitting an AGN template (broad-line 
AGNs, QSOs) are excluded since their continua will be af- 
fected by the light from the active nucleus, and therefore can- 
not be fitted with our models. In terms of area, our sample 
lies in the overlap of CFHTLS and GALEX regions (dark gray 



region in Figure [U, which contains 5878 DEEP2 galaxies. 
Other data cover this region fully, so their f ootprints are no t 
relevant. Using a technique similar to that of lBlantonl (l2006h 
we first estimate the area of the full DEEP2 EGS (Ught and 
dark gray region in Figure [1]). Our sample contains 53% of 
the total number of sources in full area, from which we ar- 
rive at an estimate of 0.31 sq. deg. for the overlap area (dark 
gray region in Figure[T]. We remind the reader that detections 
are not required in all bands as long as the object comes from 
the overlap area, so our optical sample used in SED fitting 
remains only R-bwd selected. 

We estimate galaxy parameters such as the star formation 
rate, dust attenuation, stellar mass, age, rest-frame colors and 
magnitudes, using t he ste llar population synthesis models of 
Bruzu al & Chariot! ( l2003h.^' The m ethodology is basically 
identical to that used in lSalim et al. (2007, S07), and we re- 
fer the reader to that paper for details of stellar population 
and attenuation models. Model libraries are built by consider- 
ing a wide range of star formation histories (exponentially de- 
clining continuous SF with random stochastic bursts superim- 
posed), with a range of metallicities (exact ranges are given in 
S07). Each model is dust-attenuated t o some degree accord - 
ing to a two-component prescription of^ Charlot & Falll (120001) . 
This model assumes that young populations (< 10 Myr) he 
within dense birth clouds and experience total optical depth 
of Ty. When these clouds disperse, the remaining attenuation 
is only due to the general ISM, having optical depth of /iTy, 
where /i is typically ^ 0.3. In both cases the extinction law 
of a single population is oc A""^. In our models, we allow 
for a range of tauy and /i values as described in S07. A fea- 
ture in the SED that has the greatest weight in constraining 
the dust attenuation is the UV slope, which is st eeper (the UV 
color is redder) when dust attenuation is higher dCalzetti et alj 
119941) . However, there is a significant scatter between the UV 
slope and the dust attenuation due to the differences in the SF 
history (Kong et al. 2004), which in our model is constrained 
by the inclusion of optical data. Finally, we incl ude redden- 
ing du e to the intergalactic medium, according to .Madau et alj 
(Il996l) . SFRs and stellar masses are determined assuming a 
Chabrier IMF 

The only difference in model Ubraries with respect to S07 is 
that we now construct them in < z\ih < 1 -6 range, at 0. 1 inter- 
vals in ziib, whereas in S07 libraries extended out to ziib = 0.25 
at 0.05 intervals. While library redshift resolution is finite, 
note that we use exact galaxy redshifts to scale mass and SFR 
from normalized model quantities to full absolute values. We 
test the effects of library redshift coarseness on the derivation 
of SFR and mass. On average, redshift and z\ih differ by 1/4. 
We produced a test run where we increase this difference to 
3/4 by assigning the next or the preceding library (e.g., galaxy 
at z = 0.97 is fitted with zm, = 1-1 models instead of ziib = 1.0). 
As expected, the average values of SFR and stellar mass do 
not change, but the average absolute difference is 0. 13 dex for 
SFR and 0.09 dex for stellar mass. From this we can extrapo- 
late that when redshift and ziib differ by 1/4 this deviation will 
be 0.04 and 0.03 dex, respectively. In our analysis this will 
be reflected as the small addition to the random errors. Fi- 
nally, since only models with formation age shorter than the 
age of the universe at z\ih are allowed, the number of model 

^' An update of IBruzual & Charloll 420031) models is bei ng developed to 
address issues concerning the treatment of TP-AGB stars I Maraston et all 
2006; Bruzual 2007). However, these changes will have almost no effect on 
SFRs. The systematic effect on the stellar masses will also be limited since 
we do not use the IRAC bands. 
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galaxies decreases from 10^ at zub = to 3 x 10"^ at znb = 1-6. 
Even at the high-redshift end the number of models is suffi- 
ciently large not to introduce biases in the derived parameters 
(ISalim et al.ll2007h . 

Our SED fitting involves up to 9 flux points (FUV, NUV, 
"mmt' (m*^''"'''2')cfhtls, K), their photometric errors, and the 
redshift. Photometry for various bands has been derived in 
a heterogeneous manner, but it should reflect the total fluxes 
in most bands. This will have a negligible effect on the re- 
sults. The SED fitting has one degree of freedom (scaling be- 
tween the observed and the model flux zero points). For each 
galaxy the observed flux points are compared to model flux 
points, and the goodness of the fit (x^) determines the proba- 
bility weight for the given model, and thus of the associated 
model parameters in the final probability distribution func- 
tion (PDF) of each parameter (such as the SFR, stellar mass, 
etc.). We then use the average of the probability distribution 
as our nominal estimate of a galaxy parameter and consider 
the width of the probability distribution function as an esti- 
mate of parameter error and its confidence range. In cases 
where no detection is present in a given band, that band does 
not contribute to x^. The Bayesian SED fitting performed 
here has many advantages with respect to more traditional 
maximum likelihood method. The parameter PDFs allow us 
to determine how well a given parameter can be determined 
taking into account not only the observational errors, but also 
the degeneracies among the models. For example, suppose 
that the dust attenuation and the metallicity were completely 
degenerate, i.e., that various combinations of the two produce 
identical SEDs. While the maximum likelihood will pick one 
(basically arbitrary) SED and its parameters as the best fitting, 
the Bayesian fitting will produce a wide flat PDF suggesting 
that many different values are equally probable. Similarly, the 
lack of observational constraints will also be reflected in the 
increased width of PDFs of those parameters that rely on these 
observations. While all flux points contribute to all galaxy pa- 
rameters, it is to be expected that the UV fluxes will be more 
critical in obtaining current SFRs and dust extinctions, while 
flux points red-ward of 4000A will contribute more to the de- 
termination of the stellar mass. Also, we note that despite 
the fact that our input (observed SED) contains limited in- 
formation content, one could in principle derive an unlimited 
number of galaxy parameters, since the PDF of each param- 
eter will correctly marginalize over observational and model 
uncertainties. Most of these parameters will, of course, not 
be independent, which one can again establish using (multidi- 
mensional) PDFs. Reader is refer red to S07 for further details 
about the SED fitting procedure. I Walcher et"an (l2008b . who 
use very similar model libraries and the fitting technique, also 
provide extensive discussion on the method and its robustness 
(their §2). In §|4]we will discuss in more detail the errors in 
the SFR. 

Before performing the SED fitting, we first examined color- 
color diagrams where we plot observed colors in some red- 
shift interval together with model colors corresponding to 
that redshift. We were prompted to perform these tests af- 
ter learning t hat there could be a disc repancy between the 
observed and lBruzual & Chariot! (l2003h colors in the VVDS 
sample, in the sense that models were underestimating the 
flux in the 3300-4000 A range (Wa lcher et aD 2008). By visu- 
ally comparing the locus of observed and model colors for 
various combinations of colors and redshift bins, we were 
mostly able to confirm this effect. We find the level of dis- 



crepancy (up to 0.2 mag) to be similar for both blue and red 
galaxies, which makes it unlikely to be the result of con- 
tamination fro m [OII1A3727 em ission line (emission lines are 
not included in lBruzual & Char iot 2003 modeling ) but , rather 
caused by differences in the continuum. Walcher et al. (2008J) 
use an updated version of Bruzual & Ch ariot (2003) mod- 
els (which include new prescriptions for TP-AGB stars) and 
still encounter the discrepancy. However, both the original 
iBruzual & Charlo^ (|2 003) models used here and the updated 
version used by' Walcher et al. (2008) are based on same stel- 
lar libraries which have a transition from empirical STELIB 
spectra to synthetic BaSeL spectra at 3200 A. It is beyond 
the scope of this work to try to understand the origin of this 
problem. Since at a given redshift this discrepancy would 
affect only one of our flux points, we decide to exclude 
that flux point from the SED fitting, i.e., we exclude g at 
0.3 < ziib < 0.5, r at 0.7 < zm < 0.9, / at 1.0 < ziib < 1.3, 
and z at znb = 1 -4. 

We require a minimum of three flux points for the SED fit- 
ting, though most galaxies have many more. In 336 cases this 
criterion is not met (mostly because CFHTLS magnitudes are 
not measured in spite of the fact that the object is listed), and 
we exclude these objects from further analysis. Additionally, 
we eliminate 197 objects with poor SED fits (i.e., high x^) 
whose galaxy parameters are unreliable. In S07 we discuss 
galaxies with bad SED fits and conclude that they mostly re- 
sult from bad data rather than the limited parameter space of 
the models. Thus we arrive at the final optical sample of 5345 
objects for which we obtain galaxy parameters from the SED 
fitting. 

Two parameters derived from the UV/optical SED fitting 
will feature most prominently in this work: the "current" 
SFR (i.e., the SFR averaged over the last 10*^ yr, the shortest 
timescale that can be reliably probed with stellar continuum) 
and the "age-averaged" SFR (i.e., the SFR averaged over the 
dominant population age, which depends on a galaxy and typ- 
ically ranges between 1-3 Gyr). Both will be discussed in 
more detail later Here we wish to assess the typical errors as- 
sociated with these parameters. Both SFRs are dust-corrected 
and their probability distribution functions will automatically 
reflect various sources of uncertainty. In the case of the cur- 
rent SFR the error will be dominated by the uncertainties in 
the dust correction, which we confirm by finding a strong cor- 
relation between SFR error and dust correction error. In Fig- 
ure|2](upper panel) we show the error in current SFR (fsp = 10** 
yr) as a function of rest-frame NUV-/? color, which we will 
use to select actively star-forming (blue) and quiescent (red) 
galaxies. The error equals 1/4 of the 95% confidence range of 
a PDF, which in the case of a Gaussian distribution would cor- 
respond to a standard deviation (la). The majority of galax- 
ies have errors below 0.2 dex (60%). As expected, the er- 
rors increase as one moves towards redder, less actively star- 
forming galaxies. Some galaxies, regardless of color, have an 
error in excess of a factor of 3 (0.5 dex). We find that these 
galaxies are very faint in the UV, with rest-frame FUV mag- 
nitudes fainter than 24.5. Figure |2] (lower panel) shows errors 
in age-averaged SFR. Unlike the current SFR which is mostly 
constrained by rest-frame UV, the age-averaged SFR is deter- 
mined by optical light of stars having ages 1-3 Gyr It is also 
typically less than 0.2 dex, but it is on average higher for blue 
galaxies where the optical light is fainter Unlike the current 
SFR, here the errors stay below 0.5 dex, owing to the lack of 
very faint optical sources (the sample being /?-band selected). 
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Fig. 2. — Errors in the current (upper panel) and age-averaged (lower 
panel) star formation rates, both connected for dust attenuation and plotted 
against the rest-frame color. En'ors are estimated from the width of the prob- 
ability distribution function and take into account measurement and model 
uncertainties, such as the uncertainty in the dust correction. 
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Fig. 3. — The effect of Ha contamination on the SFRs derived from the 
SED fitting. In the displayed redshift range the Ha passes through z filter (its 
relative contribution to z. is shown as a solid curve with arbitrary amplitude) 
and could thus affect the SED fitting. Shown is the difference in age-averaged 
SFRs (SFR with subscript 'a') when the z band is excluded from the fitting, 
compared to nominal fitting. Mean residuals (solid line slightly above the 
y = line) are below 0.02 dex and are not correlated with the expected Ha 
contribution. We plot only blue, star-forming galaxies for which the contri- 
bution of Ha should be the largest. 



While the stellar mass doesn't figure prominently by itself in 
this work, let us mention that the typical stellar mass errors 
are below 0. 1 dex. 
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Fig. 4. — Detection fraction of MIPS 24 ^m observations in bins of redshift 
and apparent magnitude. Gray pixels represent the detection fraction, with 
black being 100% and white being zero, except in the upper right corner 
where there are no objects. 

Galaxy SED models used in this work come from stellar 
population synthesis alone, without the inclusion of gas emis- 
sion lines. This could potentially lead to systematic effects in 
the parameters derived from the SED fitting, since the emis- 
sion lines could "contaminate" the broad-band fluxes used in 
the fitting. Typically, the most luminous emission line in our 
sample is Ha, followed by [OII]A3727. Ha becomes red- 
shifted beyond the reddest optical band (z band) at redshifts 
above 0.4, so it does not affect most of the galaxies in our 
sample. For those at lower redshift we estimate the effect of 
Ha contribution by running the SED fitting without the z band 
and comparing the resulting SFRs to those from the nominal 
fitting. The residuals of age-averaged SFRs are shown against 
the redshift in Figure [3] for blue (mostly star-forming) galax- 
ies. We present age-averaged SFRs since they should be more 
affected by the z-band flux than the current SFRs. Typical av- 
erage residuals are below 0.02 dex, and there is no obvious 
correlation with the expected relative contribution of Ha to 
z-band (solid curve, shown with arbitrary amplitude). Resid- 
uals tend to be positive, which actually corresponds to SFRs 
from the fit without the z band being larger than the nominal 
ones, the opposite from what is expected if Ha raises the z- 
band flux. As for [Oil] line, one cannot evaluate its potential 
effect on broad-band fluxes because of the unrelated issues 
with models in the 3300-4000 A range (discussed previously 
in this section). Since we already exclude from SED fitting 
the bands that sample this wavelength range, any effects of 
[Oil] wiU be removed from our nominal results. 

4. INFRARED AND OPTICAL PROPERTIES OF 24 fj.m SAMPLE 

Out of 5345 objects in the optical sample, we have 24 /im 
detections for 2430 (45%). The 24 fim imaging covers 99.5% 
of the area of the optical sample, with exposure times varying 
across the field from ^ to 19 ks (average exposure is 10 ks). 
The optical source detection efficiency grows linearly with the 
logarithm of the exposure time; it is 26% at 1 ks, and 56% at 
19 ks. In Figure|4]we plot the 24 /im detection efficiency as a 
function of redshift and apparent R magnitude. The gray scale 
is proportional to the detection fraction, with black represent- 
ing 100% and white being zero, except at bright magnitudes 
and high redshifts where there are no objects. At each red- 
shift the detection fraction increases with optical brightness, 
but for a given R magnitude the efficiency increases with red- 
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Fig. 5. — Rest-frame color distribution of the sample. Bold histogram 
shows the distribution of our optical (R < 24.1) sample (for which the SED 
fitting is peii'ormed), and the thin histogram shows the distribution of sample 
sources which are detected at 24 /tm. The dashed line is the ratio of the 
two histograms, i.e., 24 fj,m detection fraction. The 24 fim detection fraction 
peaks at the red end of the blue sequence and in the "green valley". In this and 
subsequent figures, superscript zero in NUV— i? designates rest-frame (note 
that NUV— R is not dust corrected). 



shift, which is the consequence of the detection efficiency ris- 
ing with absolute magnitude. In Figure|5]we plot the distribu- 
tion of rest-frame NllV -R colors determined from the SED 
fitting. The optical sample (bold histogram) is dominated by 
galaxies lying in the blue sequence ("(NUV-/?) < 3.5)^^. The 
peak of the red sequence (^{NW -R) > 4.5) is less obvious 
because the /?-band selection eliminates fainter red galaxies at 
higher redshifts. The thin-line histogram shows galaxies with 
24 fim detection. Again, most detections are of blue galaxies. 
The ratio of the two histograms represents the 24 ^um detec- 
tion fraction and is plotted as the dashed curve with a corre- 
sponding axis on the right-hand side. The 24 /im detection 
efficiency strongly peaks at intermediate colors, including the 
so called "green valley" (3.5 <° (NW-R) < 4.5), where it 
reaches > 80%. A sim ilar result was recently obtained by 
ICowie & Barged (l2008l) . 

Motivation for the above division into blue and red se- 
quence galaxies comes from a marked bimoda Uty in optical 
colors of local (z < 0.2) galaxies (Strateva et aLl l200ll) . where 
blue sequence galaxies have active star formation, while red 
have generally ceased forming stars. The introduction of 
UV-to-optical colors by GALEX led to a recognition of a re- 
gion in between the blue and the red sequences that was not 
prominent in optical colors CWyder et al.l l2007h . Galaxies 
that occupy this region, the green valley, acquire intermedi- 
ate colors either because they have an intermediate SF his- 
tory (transitional galaxies), or because their colors have been 
redden ed by dust and would otherwise be blue (dusty star- 
bursts) dMartin et al.l2007l;[Sarim et al.l2067b . One can distin- 
guish between the two by plotting the specific SFR (SFR/M*) 
against the color This is shown in Figure |6] where the rest- 
frame color, dust-corrected SFR and the stellar mass come 
from our UV/optical SED fitting, and the SFR is what we call 
current, i.e., averaged over 10** yr If dust reddening is mod- 
erate, there should be a correlation between the specific SFR 
(basically a ratio of recent to past S F) and the rest-frame UV- 
to-optical color dSalim et"ani2005h . One can see that this is 

Throughout the text, superscript zero designates rest-frame and not that 
the color is dust-corrected. 



\ -10 

o 

K 

M 

-2 -12 



13 



-14 




blue 



reen red 



"(NUV-i?) 

Fig. 6. — Relation between the rest-frame UV-to-optical color and the 
specific SFR of the part of the optical sample detected at 24 ^m. The two 
quantities are roughly equivalent, except that the specific SFR, derived in 
UV/optical SED fitting, is dust-corrected. Also, SFR is averaged over 10** 
yr, which we call "current" SFR (subscript 'c'). In cases of moderate dust 
UV-to-optical color correlates with the SF history of a galaxy (the trend of 
the majority of galaxies in this sample). Galaxies that scatter away from that 
relation are dusty starbursts (objects above the dashed line). Note that many 
genuine red sequence (and not just reddened), i.e. quiescent, galaxies are 
still detected at 24 /im. The classification of colors into blue, green, and 
red, and the SF histories into active star-forming, transitional and quiescent is 
based on local {z < 0.2) studies that employ photometry and spectroscopy to 
distinguish between these categories (Wvder et alJ200'i : ISalim et aU2007r) . 

the case for the majority of galaxies, especially in blue and 
red regions. However, if a galaxy has dusty star formation, 
it will have redder colors for a given specific SFR (because 
the SFR is corrected for effects of dust, while color is not). 
Such galaxies scatter above the main trend in Figure|6](galax- 
ies above the dashed line). In terms of colors, dusty starbursts 
are present in the blue sequence and the green valley, with 
the relative number of dusty to non-dusty systems peaking in 
the green valley. Here we note that even after accounting for 
dust attenuation there still exist 24 /im detections among fairly 
quiescent galaxies (log(SFR/M,j) < -11). The source of their 
mid-IR emission will be discussed in ^6] 

To infer the total infrared luminosity, LiR(8-1000/im), we 
fit 24 fj,m flux densities and redshifts to infrared SED tem- 
plates of Dale & Helou (2002).^^ These templates were nor- 
malized to follow the local //M5-cahbrated far-IR color vs. 
luminosity relation of iMarcillac et al.l (l2006l) . The assump- 
tion is that for most galaxies in our sample the mid-IR flux 
is representative of the total IR luminosity, and that one 
can use the luminosity-dependent SED models to constrain 
it. This is certainly an oversimplification. The rest-frame 
wavelength range probed by our 24/im observations (10- 
18/xm) contains many PAH lines whose relation to the mid- 
IR continuum and to the total IR luminosity may vary sig- 
nificantly compare d to th e fixed ratios assumed in the tem- 
plates aSmith et all l2007i) . Also, translating mid-IR lumi- 
nosities to total ones introduces potentially large uncertain- 
ties in the K correction. However, the use of total vs. mid- 
IR continuum luminosities is not critical in this work, and, 
as we will show in ^7.31 the results do not change if we 
use monochromatic mid-IR luminosity instead. The use of 
total IR luminosity is motivated by the commonality with 

In most of the paper total infrared luminosity will simply be called "in- 
frared luminosity", or Lir. 
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Fig. 7. — Comparison of IR luminosities obtained with various IR tem- 
plates and the 24 fim flux point. Comparison is against Ljr obtained with 
Igale & Helou (2002) templates (DH02). The difference is relatively small 
with respect to Charv & Elbaz 1 2001) templates (CEOl), but more significant 
with respect to Rieke et al. ( 2009) templates (R09), especially at higher lumi- 
nosities. We use IR luminosities based on lDale & Heloii j2002) in the rest of 
the paper 

which this measure is interpreted as a star formation rate, 
especially in high-redshift studies. Once we obtain the in- 
terpolated IR template, we calculate the total infrared lumi- 
nosity according to the relation of Sanders & Mirabel ( 1996) 
(directly integrating the SED produces very similar results, 
and Sanders & Mirabel defi nition is used as a conv ention). In 
addition to Ljr derived from lDale & Heloul (|2002[) templates, 
we additionally calculate Ljr based on luminosity-dependent 
templa tes of Chary & Elbaz (2001) and recent templates of 
iRieke et alj (l2009i). The comp arison of the two with respect 
to LiR from iDale & Heloul ( I200Z) is s hown in Figure U] For 
our sample the IR lumin osities from ICharv & ElbazI 1200 ll) 
and lDale & Heloul ( |2002|) stay within 0.2 dex of one another 
(average difference is -0.03 dex and the standard deviation 
of the ratio is 0.09 dex). At logLiR(DH02) > 11 the scat- 
ter in the ratio is very small (0.02 d ex), and the differenc e is 
almost constantly around -0.06 dex dPale & Heloul (|2002|) es- 
timate being higher). Differences with respect to lfeieke et alj 
(|2009) IR luminosities are much hig her, especially for L IRGs 
and ULIRGs (logLiR(DH02) > 11). iRieke et all (|2009|) esti- 
mates get increasingly discrepant as the luminosity increases 
(up to an order of magnitude), to the e xtent that 16% o f wha t 
are classified as LIRGs accor ding to 'Dale & Heloul (l2002h 
become ULIRGs with Rieke et al. ( 2009) templates, while 
the number of ULIRGs changes from 21 to 156. Using 
iRieke et al] (12009.) templates could possibly affect some of 
the results in our work, but also those of many other stud- 
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Fig. 8. — Infrared luminosity Ljr, derived from 24 fim flux, as a function 
of redshift. Two dashed lines define LIRGs (1 1 < logLjR < 12), and ULIRGs 
(logLiR > 12). Normal star forming galaxies (logLiR ~ 10) are detected to 
z ~ 0.7. There are very few ULIRGs in our sample. 

ies. On the other hand, the diffe rences between the other tw o 
templates are smaller (see also iPerez-G onzalez et al."2008l). 
and we adopt IR luminosities based on|Pale &Helou (200^ 
templates as our nominal values. All of these templates are 
based on local, star-forming galaxies, so they may not be en- 
tirely appropriate for high redshift galaxies such as those in 
our sample, or to more quiescent galaxies. In our analysis we 
will therefore use caution when interpreting the IR luminosi- 
ties. 

In Figure |8] Lir is shown as a function of redshift. Two 
dashed lines show regions that define LIRGs and ULIRGs. 
LIRGs start to dominate raw counts at z ^ 0.8 and we re- 
main sensitive to LIRG luminosities out to the upper redshift 
limit. We are also sensitive to normal star forming galaxies 
(Lir - 1O"'L0) to z - 0.7. The number of ULIRGs is small 
even at the highest redshi fts. This is similar to th e luminosity 
distribution presented in Le Floc'h etaP (l2005h for CDF-S. 
Since we study only 24 ^im detections with available spectro- 
scopic redshifts, we check if the redshift selection introduces 
any biases at the high-luminosity end of Lir. For this pur- 
pose we consult a cat alog of photometr ic redshifts based on 
CFHTLS photometry dllbert et alj|2006l) . and match it to op- 
tical counterparts of 24 /im sources. We then compute Lir 
based on photometric redshift. For Lir > IO'^ '^Lq the distri- 
bution of Lir of the photometric redshift sample matches the 
shape of the distribution of Lir in our spectroscopic redshift 
sample, implying no bias of the latter at high IR luminosities. 

5. INFRARED LUMINOSITY AND UV/OPTICAL SFR IN BLUE 
SEQUENCE GALAXIES 

We now have on one hand dust-corrected SFRs constrained 
from UV/optical SED fitting, and on the other hand infrared 
luminosities from 24 iim flux. We will often refer to dust- 
corrected UV/optical SFRs as SED SFRs, or just SFRs. We 
emphasize that SED fitting allows us to construct SFRs on 
various timescales, i.e., SFR averaged over some time inter- 
val. They are chosen to be fsF = 10^ 10^ 10*^ and 2 x 10*^ 
yr (averaging interval fsF ends with the epoch of the observa- 
tion, i.e., it is not centered on it). In addition to these fixed 
timescales, we also determine SFR averaged over the age of 
the galaxy, which is calculated as the total stellar mass (cur- 
rent mass plus the recycled mass as estimated in our mod- 
els) divided by the time since galaxy formation (from mod- 
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Fig. 9. — Comparison of infrared luminosities and dust-con'ected 
UV/optical star formation rates for all galaxies. Throughout the paper the 
SFRs are given in M0yr"' and for Chabrier IMF. UV/optical SFR is averaged 
over 10* yr ("current" SFR, subscript 'c'), the characteristic UV timescale. 
Numbers in lower right comer show the dispersion (cr) and the mean offset 
(A) (in dex) with respect to the 1:1 con'espondence between the SFR and 
infrared luminosity (solid line) that assumes the Kennicutt 1 1 998) conversion 
(converted to Chabrier IMF). The Kennicutt conversion is derived for dusty 
star-forming galaxies with constant SFRs over lO'-lO** yr. 

els). Because in our models the galaxy formation ages have 
a uniform distribution, i.e., they are not restricted to some 
high-redshift galaxy formation epoch, the derived formation 
age will be largely driven by the age of the dominant popula- 
tion in terms of light production (for example, blue galaxies 
will be assigned young "formation" ages regardless of their 
"real" age). We verify that there is a very tight correlation 
between the derived "formation" age and the light-weighted 
age. Therefore, what is actually measured by age-averaged 
SFR is the average SFR over the age of the dominant popu- 
lation, which for blue-sequence galaxies in our sample varies 
between 0. 1 and 3 Gyr Errors in UV/optical SFRs were dis- 
cussed in ^for the full optical sample and those conclusions 
are applicable here for the subset detected at 24 /im. 

Since the IR luminosity is usually considered in the context 
of (current) star formation rate, we begin our analysis using 
the concept of star formation rate, but extending it to include 
SFRs averaged over longer time periods. The temporal as- 
pect is not essential here. The timescales used in averaging 
the SFR simply correspond to the light emitted today by stel- 
lar populations of different ages. Therefore, what constrains 
the SFRs averaged over progressively longer timescales is the 
rest-frame luminosity at increasingly redder wavelengths. We 
will return to this relation between SFR and luminosity later. 

The timescale corresponding to lifetime of stars produc- 
ing the majority of non-ionizing UV radiation is < 10^ yr 
(lKennicutllll998h . Since young stars are typically assumed to 
dominate the dust heating at mid-IR (©, we begin by com- 
paring the infrared luminosity with SED SFR averaged over 
fsF =10** yr, which can be regarded a "current" SFR. In Figure 
|9]we plot all 2430 galaxies from our optical sample detected 
at 24 /im. The line represents the 1 : 1 correspondence between 
the SFR and the infrared luminosity assuming the Kennicutt 
conversion (converted to Chabrier IMF by applying a factor of 
1.58, S07). It is important to recall that the Kennicutt conver- 
sion applies to optically thick dusty starbursts with constant 
SF histories over 10^-10^ yr and solar metallicities. The con- 



version factor is not empirical, but is derived from population 
synthesis models. While it is not strictly appropriate to use 
this conversion for other types of galaxies (less dusty or less 
active), such practice is often encountered. This canbe some- 
what justified because, as shown in modeling of'Inoue"(2002"), 
the fortuitous cancellation of smaller dust opacity and the in- 
creased IR cirrus causes the Kennicutt conversion to also hold 
for less bursty (more normal) SF galaxies. In any case, con- 
clusions in our work are independent of the validity of con- 
versions of IR luminosity to SFR, and instead we deal with 
IR luminosities directly. In Figure |9] one sees a good over- 
all agreement between the IR luminosity and the UV/optical 
SFR. The points on average lie 0.02 dex from the 1:1 line. 
The standard deviation in the logarithm of Ljr to SFR ratio, 
i.e., the scatter around the 1:1 line, is 0.49 dex (a factor of 
three). It is the scatter that we will use as an indicator of the 
level of agreement between the IR luminosity and the SED- 
derived UV/optical SFR. 

There are a number of potential causes for the level of scat- 
ter seen in Figure |9] First, there are measurement errors. On 
average, the error in "current" SFR is 0.19 dex. The aver- 
age error in 24 /^m flux measurement is 6% (0.02 dex), which 
is negligible in comparison. We expect some error from the 
extrapolation of the observed 24 /im (r est frame 10-18 jum) 
luminosity to total IR luminosity (e.g, Le Floc'h et alj|2005l 
claim a factor of three error). As discussed, the two frequently 
used sets of IR SED templates already produce a scatter of 0. 1 
dex in their estimate of Lir. Another source of scatter could be 
from from the inclusion of all galaxies in our sample, includ- 
ing many red galaxies with older stellar populations and not 
much current star formation. In order to compare UV/optical 
SFR and Ljr with an assumption that IR arises from star for- 
mation, one needs to limit the sample to actively star forming 
galaxies. Following discussion in ^it would be appropriate 
to base such selection on a dust-corrected quantity such as the 
specific SFR. However, it is more intuitive to use rest-frame 
color instead. Taking blue galaxies ("(NUV-/?) < 3.5) will 
select most actively star-forming galaxies (including a large 
number of dusty starbursts), while not allowing galaxies with 
more quiescent SF histories (Figure|6l). In Figure[TOl we now 
compare Ljr and SED SFR of blue galaxies alone. The SFR 
averaging timescale is still 10** yr The scatter in Lir to SFR 
ratio is 0.42 dex, or 16% smaller than in the full sample. The 
scatter was reduced by the removal of red galaxies. This re- 
duction cannot be attributed to slightly smaller SED SFR er- 
rors: 0.17 dex for blue galaxies vs. 0.19 dex for the full sam- 
ple. 

Throughout the LIRG range of luminosities the agreement 
between the IR luminosity and UV/optical SFRs is relatively 
good, albeit with large scatter. This implies that the LIRGs 
in the redshift range studied here cannot be optically thick at 
UV and optical wavelengths, thus allowing us to use the stel- 
lar continuum to deduce the SFR and other parameters, such 
as the stellar mass. This, of course, depends on our ability 
to obtain reliable rest-frame luminosities and dust attenuation 
estimates. One also sees that the slope betwe en IR luminos- 
ity and UV/optical SFRs is steeper than the Kennicutt (1 19981) 
relation. This is fully expected. The Kennicutt relation ap- 
plies to galaxies in which a large fraction of stellar emission 
is absorbed by dust. This will be less the case for galax- 
ies with smaller SFRs, w hich have smaller dust attenuations 
(Wang & Heckman 1996!). 

Next we explore how the overall scatter in Lir to UV/optical 
SFR ratio changes if SFR is averaged over timescales other 
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Fig. 10. — Comparison of infrared luminosities and dust-corrected 
UV/optical star formation rates for blue-sequence galaxies. Timescale for 
SF is still 10^ yr, and the Kennicii^ fT993) conversion (derived for dusty star- 
forming galaxies with constant SFRs over lO'-lO* yr) is shown as the solid 
hne. Removal of red galaxies reduces the scatter in correlation. Numbers 
have the same meaning as in Fig.|9] 
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Fig. 11. — Comparison of infrared luminosities and dust-corrected 
UV/optical star formation rates for blue-sequence galaxies, where SFR is now 
averaged over the galaxy age (i.e., the age of the dominant population, 0.1-3 
Gyr old, subscript 'a'). The correlation between Ljr and SFR is better (the 
scatter is smaller) than in Figure [Tol where SFR was averaged over 10** yr. 
The en'or in the SFR determination is similar here as it is for SFR averaged 
over 10* yr (~ 0.2 dex), and does not dominate t he scatter with respect to 
LiR. There is a departure with respect tol KennicuE (,1998) conversion which 
is derived for dusty star-forming galaxies with constant SFRs over 
yr, because most galaxies have declining SF histories so the SFRs averaged 
over longer timescales are on average higher than the SFRs averaged over 
10** yr. Numbers have the same meaning as in Fig.|9] 

than fsF = 10^, still restricting our focus to blue galaxies. We 
again emphasize that averaging SFR over shorter or longer 
time periods is a way to probe the connection of today's stars 
having different range of ages using the SFR concept, and 
does not imply that the past episodes of star formation di- 
rectly affect the IR luminosity that we see today. We start 
from fsF = 10^ yr, which can be considered an "instantaneous" 
SFR, and find scatter to be 0.43 dex, slightly larger than for 
fsF = 10^. Note that UV/optical SED fitting does not constrain 
the SFR averaged over such short timescales very well, and 




LSF residuals 

Fig. 12. — Distribution of the residuals of the linear fit of infrared luminosi- 
ties and dust-corrected UV/optical star formation rates. Upper panel shows 
residuals with respect to SFRs averaged over 10** yr, and the lower with re- 
spect to SFRs averaged over the population age. A Gaussian is fit to each 
distribution (dashed curve), and its width (ctghuss) is displayed. The Gaus- 
sian of the residuals with respect to SFRs averaged over the population age is 
narrower The horizontal axis is in decades (dex). 

the increase in scatter is simply the result of the poorer qual- 
ity of the SFR measure over this timescale. But for the two 
longer timescales, fsF = 10^ yr and fsF = 2 x 10^ yr, the scat- 
ter decreases, to 0.39 and 0.37 dex, respectively. We obtain 
yet smaller scatter, 0.34 dex, when we consider SFR averaged 
over the age of the dominant population in the galaxy, shown 
in Figure[TT] This represents a 20% reduction in Ljr vs. SFR 
scatter compared to fsF = 10^ SF timescale. The values of 
scatter we give here are for the Lir to SED SFR ratio, but very 
similar answers are obtained if we consider a scatter around 
the best linear fit. 

One may get an impression that most of the reduction in 
scatter as we go to longer SFR-averaging timescales is due to 
fewer outliers. To check this, in Figure [12] we fit Gaussian 
functions to residuals around the linear fits, and find that the 
width of the Gaussian (which is not dominated by outliers) 
decreases similarly as the overall scatter, indicating that the 
decrease in scatter is not due to the decrease in the number of 
outliers. 

The most pressing concern with the above result is that the 
reduction in scatter when comparing Lir to SFR over progres- 
sively longer timescales is an artifact of the SED fitting pro- 
cedure, such that it simply reflects the precision with which 
we are able measure SFRs (and dust corrections) at different 
timescales (or alternatively, wavelengths)? The average for- 
mal error in our SFR measurements is between 0.14 and 0.18 
dex for SFRs averaged over timescales 10^ yr and longer The 
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small value of the error and its small variation for different 
timescales implies that the SFR uncertainties are not modulat- 
ing the level of correlation with Lir, i.e., the change in scatter 
is not driven by errors in the measurement of the SFR. While 
the above is true on average, in Owe saw that in some cases 
the error in SFR averaged over 10** yr can get relatively high. 
Removing all galaxies with error larger than 0.2 dex leads to 
some reduction in scatter with respect to Lir, but is still larger 
than the scatter between Ljr and age-averaged SFR. The same 
is true if we limit the sample only to objects where the error 
in SFR averaged over 10'^ yr is smaller than the error of age- 
averaged SFR. 

To further test if we are able to reliably measure the SFR on 
10** yr timescale, we run simulations described in Appendix 
A. From those we conclude that if the IR luminosity were in- 
deed the reflection of the current SF, then our UV/optical SFR 
averaged over 10** yr would measure it with a smaller scat- 
ter than an UV/optical SFR averaged over any other longer 
timescale. 

Finally, we notice that the average offset with respect to 
Kennicutt conversion is larger when age-averaged SFRs are 
plotted instead of the current ones in Figure [TT] This should 
not be surprising since the Kennicutt conversion was cali- 
brated assuming current (< 10** yr) SFR, and in general the 
SFR averaged over longer time periods will be higher than the 
current one (because most galaxies have declining star forma- 
tion histories). Additionally, the correlation is now steeper 
(the slope from the bisector linear fit was 1.06 for fsF =10** 
and is now 1.20). Since in both cases one has the same 
LiR, the change in slope has to be the result of a differen- 
tial change in SFR between the two averaging timescales for 
galaxies with low and with high IR luminosities. As can be 
seen from Figure |8] galaxies with Lir < lO'^L© are detected 
only at redshifts below 0.5. Galaxies that we see at this lower 
redshift will on average be older than the galaxies observed 
at higher redshift, when the universe was younger. For the 
same rate of SF decline, galaxies that had more time to evolve 
(Lir < IQ^^Lq galaxies at z < 0.5) will show greater change 
between the age-averaged SFR and the current SFR than the 
younger galaxies (more luminous). This moves Lir < lO'^L© 
galaxies more to the right in Figure[TT]than the more luminous 
ones, producing a steeper slope. 

5.1. Lir and UV/optical SFR: dependence on galaxy color 

The source of the IR luminosity will generally not be the 
same in galaxies with different dominant stellar populations. 
Therefore, we now explore the strength of the Lir vs. SFR 
correlation, not only as a function of timescale, but also for 
galaxies split into various color bins. We now use rest-frame 
B-R color because it somewhat better discriminates the pop- 
ulation age of star-forming galaxies than NUV-/?. In Figure 
[T3]we plot the scatter of the logarithm of Lir to SFR ratio as a 
function of timescale over which the SFR was averaged. Each 
panel displays the same relative dynamic range in the verti- 
cal axis. The top three panels show blue-sequence galaxies, 
while the bottom panel contains the green valley and the red- 
sequence galaxies. Open squares represent the scatter for cor- 
responding fixed timescales, while the filled square is the scat- 
ter in age-averaged SFR, plotted at the position of the average 
age of galaxies in that color bin. Even for the bluest galax- 
ies (^Hfi-R) < 0.3, top panel), which have a large fraction of 
recent star formation, the scatter is smallest at a timescale of 
10^ yr, rather than the 10** yr UV timescale. In subsequent 
redder color bins the best correlation with Lir is always for 
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Fig. 13. — Correlation between infrared luminosities and dust-corrected 
UV/optical star formation as a function of SF averaging timescale, for galax- 
ies having different rest-frame colors. The top three panels correspond to 
blue-sequence galaxies, and the bottom panel contains green valley and red- 
sequence galaxies. The correlation generally gets better (the scatter de- 
creases) as the timescales increase. Filled squares correspond to the age- 
averaged SFR and they are plotted at /jp corresponding to the average age of 
galaxies in a given color bin. Unlike other timescales, the SFR averaged over 
lO' yr is poorly constrained in UV/optical SED fitting and would increase the 
intrinsic scatter in Lir/SFR congelation. Each panel displays the same relative 
dynamic range in vertical axis. Photometric bands characteristic for a given 
timescale are given above the plot. 

age-averaged SFR, i.e., on timescales of ^ 2-3 Gyr. We tried 
to identify a galaxy population for which the IR luminosity 
would best match the short timescale of 10** yr. We looked 
at the Lir/SFR scatter in bins of galaxy stellar mass (M»), 
specific SFR (SFR/M,) and the age of the most recent burst. 
We find that IR best correlates with 10** yr timescale only for 
logM* < 8.5 galaxies. These are blue compact dwarfs that we 
can detect only out to z ~ 0.4, so the result is based on a small 
number of objects. 

As explained previously, we begin our analysis using the 
concept of SFRs averaged over various timescales. However, 
this was simply a convenient way to probe today 's stellar pop- 
ulations of different ages. Given that the flux that is respon- 
sible for dust heating must be produced at the present time, 
a quantity that will be more fundamentally correlated to Lir 
is some UV or optical luminosity. For every SFR-averaging 
timescale fsF there is a characteristic (rest-frame) wavelength 
at which the population with the age fsF dominates. In Figure 
[T3]we show which of the bandpasses (FUV, NU V, U, B and 
V) cor respond to various ages, extrapolated from lO'Connelll 
(1990). The FUV will be dominated by stars having ages 
< 10'^ yr, so the equivalent to "current" SFR will be the dust- 
corrected FUV luminosity. Similarly, the equivalent for SFR 
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Fig. 14. — Comparison of infrared luminosities and dust-corrected FUV (left) and B-band luminosities (right) of blue-sequence galaxies. Dust corrected FUV 
and B luminosities have been derived from the UV/optical SED fitting. Average correction is 2.2 mag in FUV and 0.80 mag in B. The scatter around the least 
square linear fit is given in the lower right corner, and is smaller for the 6 band. Also given is the corre lati on coefficient. These correlations are equivalent to 
correlations between Ljr and the SFR averaged over the short and long timescales presented in Figs. llOl and ll II but are more fundamental in the sense that they 
tie LiR with present-day source of IR heating. 



averaged over 1-3 Gyr will be optical luminosity (U, B or V), 
also corrected for dust. The bandpass that corresponds to a 
timescale with the least scatter in top panel of Figure[T3]is be- 
tween rest frames U and B, and around V band for somewhat 
redder blue-sequence galaxies (second and third panels). Now 
we again show IR luminosities, but against dust-corrected UV 
or optical luminosities instead of the SFRs. Fi gure[T4lpresents 
a comparison of Ljr and dust-corrected FUV (left) and B 
(right) luminosities, again for blue-sequence galaxies. As 
in the case of SFRs, the dust correction for FUV and B band 
luminosities come from our SED fitting, and it is on average 
2.2 mag in FUV and 0.8 mag in B. These figures are equiva- 
lent to those that showed Lir and SFR (Figs. [Tol and [TTll. and 
the arguments that applied for the robustness of SFRs (Ap- 
pendix A) apply here for FUV and B luminosities. The scatter 
is smaller against Lb^cou-. Formal scatter around the linear fit 
is 0.39 dex for FUV luminosity, vs. 0.32 for Z?-band luminos- 
ity (0.36 and 0.30 dex when 3(j outliers are excluded). Pear- 
son correlation coefficient provides another way to measure 
the strength of a correlation. It is 0.80 for FUV luminosity 
and 0.86 for B-band luminosity. Note that in order to have a 
meaningful comparison with IR luminosity, the UV or the op- 
tical luminosity needs to be appropriately corrected for dust. 
In absence of dust correction, the correlation coefficient be- 
tween FUV luminosity and Lir drops to 0.59 and between B- 
band luminosity and Lir reduces slightly to 0.84 (since dust 
corrections in B are smaller). 

The linear fit without 3a outliers for Z?-band (A = 4360 A) 
is given by 

logLiR= 1.125(0.013)logLB.eorT-l. 102(0.142), logL^.^n > 

(1) 

where all luminosities are in Lq. The fit is constructed for 
blue-sequence galaxies ("(NUV-7?) < 3.5). The values of 
parameters of the fit depend slightly (±0.05 in slope) on the 
exact color cut. The appropriate dust correction for Z?-band 

Results for B-band luminosity are very similar to those for U or V, but 
we use B since it is the most common band used for galaxy magnitudes. 



can b e obtained by scaling the attenuation in FUV using the 
mean lCharlot & FaUl (l2000l) extinction law for age of ^ 1 Gyr 



Ab = 0.37Afuv, 



(2) 



where Apuv is preferably obtained from full SED fitting, or 
alternatively using the UV slope relation given in Equation 
BlKbut see discussion in Appendix B). 

For some purp oses, one may prefer a bisector linear fit 
(llsobe et al.lll990 ). given by 



logLiR= 1.275(0.016) logLB,eon-2.668(0. 1673), logL^.^o 

(3) 

which was constructed from all points in Figure [14] (right). 

The above correlation (Eqn. [T] demonstrates that one can 
essentially estimate, to within a factor of two, the total IR 
luminosity from UV/optical photometry alone (i.e., the dust 
corrections are also constrained only using UV/optical SED). 
While this may not be the case for every type of galaxy at 
any redshift, it appears true for normal and IR luminous star 
forming galaxies over the z < 1 .4 redshift range studied here. 
Note that we obtain Lir from 24 /im flux, and that the corre- 
lation could perhaps be even tighter with a better estimate of 
total Lir, one that employs longer wavelength IR data and/or 
more accurate SED templates. This will be addressed in fu- 
ture work. 

5.2. Lir and UV/optical SFR: dependence on redshift 

So far we have investigated the correlation between Lir and 
, SFR for various samples irrespective of their redshift. Since 
the observed 24 /im flux corresponds to 10-18 /^jm rest-frame 
wavelength range that contains both mid-lR continuum and 
strong PAH lines, one would like to learn if there are any 
systematic differences in the Lir vs. SFR correlation at dif- 
ferent IR wavelengths, i.e., redshifts. Still focusing on blue, 
star-forming galaxies, we find that the Lir vs. SFR relation in 
different redshift bins follows the same trend as found for the 
entire sample: Lir correlates better with SFR averaged over 
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Fig. 15. — Comparison of infrared luminosities and dust-corrected B-band luminosity in different redshift bins (for blue-sequence galaxies). Dust corrected B 
luminosity has been derived from the UV/optical SED fitting. Numbers in lower right comers show the dispersion around the least square Unear fit (in dex) and 
the coiTelation coefficient. The dashed line represent the bisector linear fit for the entire sample, and is repeated from panel to panel as a guiding line. 



galaxy population age than over any shorter timescale. Equiv- 
alently, and more fundamentally, Lir correlates better with 
dust-corrected optical luminosity than with dust-corrected UV 
luminosity. This is shown in Figure [15] where we present Lir 
vs. dust-corrected Z?-band luminosity of blue galaxies, split 
into six 0.2-wide redshift bins in the 0.2 < z < 1.4 range. 
The upper number in the lower right corner of each panel is 
the standard deviation around the least square linear fit, and 
the lower number the correlation coefficient. The scatter is 
roughly the same in all redshift bins, indicating that at this 
level of precision the entire 10-18 fim wavelength range cor- 
responds equally well to the UV/optical dust-corrected lumi- 
nosity. The PAH features at 11.3 and 12.7 fim would be sam- 
pled in 0.8 < z < 1 .0 and 1 .0 < z < 1 .2 redshift bins, respec- 
tively, and there we see a slight increase in scatter compared 
to other redshift bins. In each panel we repeat the bisector 
linear fit obtained for the full sample (Eqn.|3]as a dashed line. 
From that one can see that the slope appears to get steeper at 
higher redshifts. Rather than assuming that the intrinsic lin- 
ear relation changes at different redshifts, it is possible that 
this is because the intrinsic relation is not linear Then, since 
at different redshifts one samples different ranges in luminos- 
ity, segments of a curve will appear as linear relations with 
different slopes. Also, it is plausible that the conversion from 
rest-frame mid-IR flux to to Lir has wavelength-dependent 
systematics. 

One may wonder if the improvement in the correlation be- 
tween UV/optical luminosity and Lir as we go to redder op- 
tical luminosities (Fig. [T4l i may in fact reflect a more funda- 
mental correlation of Lir with stellar mass. On the one hand 
this is not expected because dust heating should correlate with 
some form of present flux, whatever from younger or older 



stars (or some mix of two), and not on mass that includes 
stars that cannot contribute significantly to the IR. On the 
other hand, we know that for actively star forming galaxies 
the sta r formation rate and stellar mass are tig htly correlated 
(e.g., Boselh et al.ll20011:lBrinchmann et al.ll2004) . When we 
compare the Lir of blue galaxies vs. their current stellar mass, 
the scatter around the least square linear fit is 0.41 dex, and 
the correlation coefficient is 0.76, which is weaker than what 
we found when comparing Lir to either FUV or Z?-band lumi- 
nosity (corrected for dust), with correlation coefficients 0.80 
and 0.86, respectively. Very similar results are obtained when 
we substitute the current stellar mass with the estimate of to- 
tal stellar mass formed over the galaxy lifetime, i.e., the mass 
that includes recycling. However, this is not the full story. 
There is evidence that the SFR vs. mass relat ion evolves with 
redshift (iPapovich et all 120061: iNoeske et alll2007b ). so for a 
given mass galaxies at different redshifts will have different 
Lir. Indeed, if we split the sample in 0.2 wide redshift bins 
(as in Fig. [TST i, the correlation between IR luminosity and 
the stellar mass improves, and is comparable to that between 
IR luminosity and the dust-corrected B-band luminosity in a 
given redshift bin ( Fig. [TsT l. Since the mass measurement in 
the SED fitting is constrained by very similar information that 
constrain the optical dust-corrected luminosity, this similarity 
between and M^xoit should not be surprising or considered 
or fundamental, but instead reiterates the connection between 
the IR emission and the stars other than the very young ones. 

6. INFRARED LUMINOSITY AND UV/OPTICAL SFR IN GREEN 
VALLEY AND RED SEQUENCE GALAXIES 

The analysis presented so far has focused on blue sequence 
galaxies, for which it was reasonable to assume that IR emis- 
sion would be strongly related to active star formation. The 
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Fig. 16. — Comparison of infrared luminosities and dust-corrected 
UV/optical star formation rates for green valley and red sequence galaxies. 
Timescale for SED SFR is 10** yr (subscript 'c'). Galaxies classified as dusty 
stai'bursts based on Figure|6]are shown as filled squares. Other red galaxies 
mostly lie above the Kennicutt ( 1998) conversion (solid line) indicating that 
the current SF is not the primary source of IR luminosity. 

picture becomes more complex as one moves away from the 
blue sequence into the green valley and the red sequence. 
While some galaxies in the green valley will simply be red- 
dened actively star-forming galaxies, others will have such 
colors because they have little ongoing star formation (cor- 
responding to galaxies above and below the dashed line in 
Figure |6]l. We should expect older stellar populations to con- 
tribute more to the IR emission in the latter group. This will 
be even more the case for red sequence galaxies, which have 
little or no current star formation. In Figure |6] we saw that 
there exist 24 /im-detected galaxies well into the red sequence 
(as red as any galaxy in our optical sample). Their low spe- 
cific SFRs indicate that these galaxies are intrinsically very 
red and not just dust reddened. In Figure [16] we compare IR 
luminosity and the SED-derived current SFR (fsp =10^ yr) 
for red galaxies (which includes the green valley and the red 
sequence). We distinguish between dusty starbursts (galaxies 
above the dashed line in Figure|6] plotted as squares) and reg- 
ular red galaxies (dots). The two groups occupy distinct loca- 
tions. Dusty starbursts have high UV/optical SFRs: above 10 
MQyr"', and in some cases approaching 1000 MQyr"'. They 
lie close to the 1:1 Kennicutt (1998) conversion between Lir 
and SFR. This is expected if Lir in dusty starbursts is due 
to SF. Actually, galaxies with the most intense SF have some- 
what lower Lir than the expected, ULIRG levels. For such ex- 
treme cases it is possible that the SED fitting overpredicts the 
dust correction (but we cannot exclude that IR luminosities 
are perhaps underestimated). Non-dusty red galaxies (dots) 
have lower SFRs and lie above the Lir-SFR conversion. This 
means that Lir is not powered by the current SF. At each 
UV/optical SFR there is a wide range of IR luminosities. This 
again speaks of a disconnect between Lir and SFR. 

Non-dusty red galaxies are the main subject of the analy- 
sis in this section. Can we explain the presence of 24 /im 
emission and the derived Lir luminosities in these galaxies 
only with stellar emission, which by necessity (since there 
is little current SF) will mostly come from intermediate and 
old stellar populations? Do we see evidence that some other 
dust heating mechanism, such as an AGN, may be present 
in these galaxies? In Figure [17] the IR luminosity for galax- 




-8 -9 -10 -11 -12 -13 
log {SFRJM,) 



Fig . 17. — IR luminosity of galaxies with different specific SFRs. Actively 
star forming galaxies are to the left of the dashed line, while those to the 
right are transitional or quiescent and they con'espond to green valley and 
red-sequence galaxies. LlRG-like luminosities (logLjR > 11) can be found 
even among some very quiescent galaxies. 

ies with different (current) specific SFRs is shown. Since, 
unlike color, the specific SFR is corrected for dust, this plot 
enables us to place red dusty starbursts together with blue ac- 
tively SF galaxies (left, log(SFR/M*) > -10), and separate 
more quiescent, red galaxies (right, log(SFR/M)* < -10) for 
which we investigate the source of IR emission. We see that 
galaxies with LIRG-like luminosities are present well beyond 
the region of actively star-forming galaxies, with some having 
specific SFRs as low as log(SFR/M*) = -11.5, which corre- 
sponds to rest-frame color of N UV-i^ = 5, the co lor of the 
bluest nearby elliptical galaxies dDonas et alj|2007 ). One can 
be concerned that the use of IR templates based on actively SF 
galaxies to derive Lir for these more quiescent objects may 
not be appropriate. This is entirely possible. However, as 
in the previous analysis, we will assume that this (commonly 
used) procedure is correct and then draw consequences. At 
each specific SFR there is a wide range of Lir, especially for 
active galaxies. This is mostly the consequence of a wide 
range of masses probed at each specific SFR. 

In order to establish if the IR luminosities that we see in 
red galaxies can be produced by stars alone (of any age) we 
perform the following exercise. The UV/optical SED fitting 
allows us to estimate the total amount of stellar luminosity ab- 
sorbedhy the dust (Cortese et al. 2008; da Cunha et al. 2008). 
According to the dust model of lCharlot & Fall ( 2000), this en- 
ergy (dust luminosity), will come from birth clouds surround- 
ing young stars (< 10 Myr old) and from the ISM heated by 
stars of intermediate and older age. In the case when there is 
no non-stellar source of IR emission, the Lir should match, or 
at least not significantly exceed the Ljust estimate. Estimating 
the amount of dust extinction in quiescent galaxies from the 
UV/optical SED fitting will be more uncertain than in actively 
star-forming galaxies, as suggested in Figure [2] Nonetheless, 
we expect that the dust luminosity derived from the SED fit- 
ting should on average be correct, which therefore allows us to 
check the energy budget. In Figure [18] we present the ratio of 
the dust luminosity (Ldu.st) derived from the UV/optical SED 
fitting to the observed IR luminosity against the current spe- 
cific SFR. Objects to the right of log(SFR/M*) = -10 line are 
red, quiescent (or transitional) galaxies. Ratio of Ljust to Lir 
of unity means that the energy that is estimates to be absorbed 
in the UV/optical part of the spectrum equals the energy re- 
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Fig. 18. — Ratio of dust luminosity to the observed IR luminosity against 
the current specific SFR. L^ust is stellar energy absorbed by the dust and is 
derived from the UV/optical SED fitting. Most galaxies, both actively star 
forming (left of the dashed line) and quiescent (together with transitional; 
right of the dashed line) have ratios around unity (horizontal line), indicat- 
ing that the dust absorption of stellar light can on average account for the 
observed Ljr. The two thick line represent ±1 cr range of error in the ratio. 

emitted in the IR. To see what range of values are consistent 
with the ratio of 1, we calculate average 68% confidence range 
of the ratio (two thick lines) from PDF errors of Ldust, with an 
ad hoc 0.3 dex error for Lir added in quadrature. Most of the 
actively SF galaxies fall in the region consistent with the ra- 
tio of unity except those with high specific SFRs. As already 
mentioned, for these objects SED SFR (and therefore Ldust) 
may be overestimated, or their Lir underestimated. Red qui- 
escent galaxies have a larger scatter of Ldust/i'iR ratios, which 
is not surprising given the higher uncertainties in estimating 
^dust from the SED fitting for these galaxies (the thick lines). 
Again, most galaxies lie within the ±la range around unity. 
From this we conclude that the dust heated by stellar pop- 
ulations is roughly sufficient to account for the observed IR 
luminosity even for relatively quiescent galaxies with LIRG- 
like IR luminosities. Consequently, we conclude that there 
cannot be a large population of (presumably obscured) AGNs 
which would significantly raise Ljr and skew the ratio below 
unity (we will next see that AGN may be affecting Lir, but 
only at a moderate level). Given the low levels of current SF 
in transitional and quiescent galaxies, one must conclude that 
intermediate and older stellar populations produce the bulk of 
the IR emission (see also Figure |23]|. 

Regardless of the arguments laid out above, one would still 
like to test directly whether the presence of AGNs has a sig- 
nificant effect on the mid-IR emission in our sample, espe- 
cially among the more quiescent galaxies. Obtaining a full 
census of AGNs in our sample is not straightforward. First, 
our UV/optical SFRs can be derived only for galaxies where 
an AGN has no effect on the UV continuum, which is why we 
have already excluded several tens of broad-line AGNs (type 1 
AGN), as identified from the spectra. To identify narrow-line 
AGNs (typ e 2 AGN) using the BPT emission line classifica- 
tion (BPT, Baldwin et al. 1981) requires spectra that cover rest 
frame range from 4800-6600 A. For our spectra this is pos- 
sible only in a very small redshift range (0.33 < z < 0.38). 
However, coupled with information on stellar mass, some 
fraction of galaxies lying in the AGN parts of the BPT diag- 
nostic diagram can be distinguished even in single-axis pro- 
jections of the diagram, i.e., using one line ratio. Using simi- 



FlG . 19. — Ratio of dust luminosity and the observed IR luminosity against 
the current specific SFR for galaxies identified as AGNs. AGNs are identified 
from X-ray detections (star symbols) and from optical emission lines (open 
diamonds). Relatively large number of AGNs, especially optically identified, 
have transitional or quiescent SF histories (right of the dashed line). Most 
AGNs have L^ust to Lir ratio consistent with unity (the shaded region between 
two thick lines represent ± 1 cr range of eiTor in the ratio, based on the full 
sample (Fig. [18)). This suggests that AGN heating of the dust is on average 
not very significant. Exceptions may be AGNs lying below the lower thick 
line with -11 < log(SFR/A/,) < -10. 

lar criteria as lWeiner et al.l (120071) . we select AGN candidates 
at z < 0.38 by requiring the flux ratio log([NII]/Ha) > -0.2 
and steflai- mass logM* > 9.5, and at 0.34 < z < 0.82 by se- 
lecting log([OIII]/H/3) > 0.7 together with logM* > 10.2. 
These criteria select a total of 35 type 2 AGN candidates de- 
tected at 24 /im, which we call "optical AGN". Additiona lly, 
we use a catalog of Chandra sources jLaird et aTl 120091) to 
identify 74 X-ray AGN candidates detected at 24 /xm. The 
majority of X-ray sources in the EGS are believed to be AGNs 
or have an AGN component (Laird et al. 2009). We plot the 
ratio of Ldust to Lir against the current specific SFR in Figure 
[19] coding points by AGN type. First, we notice that optical 
AGNs (open diamonds) are almost exclusively transitional or 
quiescent objects. This agrees with the results of local stud- 
ies where there appears to be a relation betwee n optical AGN 
and SF quenchin g (Kauffmann et al. 2003; S alim et alj|2007i : 
[Graves et al1l2007:,Schawinski et al..,200 7). X-ray AGNs are 
additionally present among the galaxies with higher specific 
SFRs, but only up to log(SFR/M,j) = -9, which again may 
be related to their role in SF quenching. Most AGNs have 
Ldust/i'iR consistent with unity (the ±1 cr range, shaded, is re- 
peated from Figure [TSl). If AGN contributes significantly to 
Lir this would be reflected in Ldust/^iR below that of non- 
AGN galaxies. While this is generally not the case, there is 
a group of AGN at -11 < log(SFR/M,s) < -10 where in in- 
dividual cases the AGN contribution to Lir may be around 
90%. 

Next we try to estimate the fraction of Lir that is on average 
attributable to AGNs. Figure|20]displays average IR luminosi- 
ties in bins of current specific SFR for three classes of galax- 
ies: 1) non-AGNs (shaded region between two thick lines), 
2) X-ray AGNs (star symbols) and 3) optical AGNs (open di- 
amonds). Error bars give the error of the mean in each bin. 
On the actively SF side we have only X-ray AGNs, and their 
average Lir is consistent with those of non-AGNs. On the 
quiescent side X-ray AGN have Lir up to 0.2 dex higher than 
non-AGNs, although it is only at log(SFR/M*) = -10.3 that 
the excess is somewhat significant. Optical AGNs are similar 
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Fig. 20. — Average Ljr of AGNs and non-AGNs in bins of specific SFR. 
Plotted are: 1) non-AGNs (shaded region between two tliick lines), 2) X-ray 
AGNs (star symbols) and 3) optical AGNs (open diamonds). Error bars give 
the error of the mean in each bin. The error range for non-AGN is given 
by the shaded region. X-ray AGNs are somewhat more luminous than non- 
AGNs around log(SFR/M, ) = -10.3. Two sets of error bars are slightly offset 
between each other in horizontal direction for clarity. 



Fig. 21 . — Infrared luminosity distribution of type 2 AGNs selected using 
a single-line BPT diagram (optical AGN, thick histogram), and of the control 
sample of non-AGNs (thin histogram). Non-detections at 24 ^m are plotted 
at logLiR = in both case. Each group contains 47 objects. Non-AGNs were 
selected to have similar masses and specific SFRs as AGNs. There is no 
significant difference in IR luminosities between the two groups, or in the 
number of 24 fim non-detections (12 for AGN, 10 for non-AGN. 



to X-ray AGNs at transitional specific SFRs, and then sig- 
nificantly lower than non-AGNs for low specific SFRs, most 
probably because these are low-power AGNs such as LIN- 
ERs. 

The above procedure has a drawback that if AGN selection 
is biased with respect to Ljr, their average Lir will be off. 
Thus we append it with the following test. For each group of 
AGNs (optical and X-ray) we select a control group of non- 
AGNs with similar properties. This needs to be done for all 
AGN (47 optical and 86 X-ray) regardless of whether they 
have been detected at 24 /im. For each optical AGN we select 
an object from the same redshift range (z < 0.38 or 0.34 < z < 
0.82) such that the emission lines do not indicate an AGN, and 
with a matching stellar mass and specific SFR. The matching 
object is defined as the one that minimizes the "distance" D 
in the stellar mass-specific SFR space: 



= (AlogM,f + c^A(log(SFR/M,)f 



(4) 



where c is a "scaling" ratio between logM, and log(SFR/M*), 
which we nominally take to be 3 based on the range of these 
quantities in our sample of AGN. Similarity in stellar mass 
and current specific SFR will ensure similarity in^rnany other 
non-AGN characteristics as well (ISchiminovich et al.ll2007l) . 
For X-ray AGNs we select matching non-AGNs such that they 
are not detected in X-ray, have a redshift within 0.2, and mini- 
mize Equation]?] For both samples the same non-AGN match 
is allowed to appear more than once. 

In Figure ]2T] we compare the distribution of IR luminosi- 
ties for optical AGNs (thick histogram) vs. non-AGNs (thin 
histogram). Objects not detected at 24 /im are plotted with 
logLiR = 0. The two distributions are quite similar, including 
the similar number of 24 /im non-detections. The average Lir 
of the AGN is 0. 16 dex higher than that of the non-AGNs (for 
the part of the sample where AGN and non-AGN are both de- 
tected at 24 /xm). However, the average stellar mass of the 
AGN sample is also slightly higher (0.24 dex), making the 
difference in Lir less significant. Results are similar when we 
choose c = 2 or 4 in Equation ]4] From this we conclude that 
optical AGNs are drawn from the same underlying IR popula- 
tion as non-AGNs. A similar comparison is shown for X-ray 
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Fig. 22. — Infrared luminosity distribution of X-ray detected AGNs (thick 
histogram), and of the control sample of non-AGNs (thin histogram). Non- 
detections at 24 /im are plotted at logLiR = in both case. Each group con- 
tains 86 objects. Non-AGNs were selected to have similar masses and current 
specific SFRs as AGNs. Non-AGN are more frequently (24 vs. 12 galaxies) 
not detected at 24 fim, making them on average somewhat less luminous in 
IR. 

AGNs in Figure ]22] If both an AGN and a matching non- 
AGN are detected in 24 /xm, their Lir (and stellar mass) are 
on average quite similar, as was the case with optical AGNs. 
However, while 12 X-ray AGN are not detected in 24 /im, 
this number jumps to 24 for the non-AGN control group. If 
we assume that each non-detection has Lir corresponding to 
the detection limit at the given redshift, we get that the aver- 
age Lir of X-ray AGNs is 0.23 dex higher than of non-AGNs 
(both AGNs and non-AGNs have the same stellar masses). 
This result is robust if we choose c = 2 or 4 in Equation]!] 

While our sample of AGN is not large enough to draw firm 
conclusions, it appears that AGN are not a significant contrib- 
utor to mid-lR luminosities in the general case, i.e., in sam- 
ples that have an optical selection, such as ours. Where their 
presence could be detected, especially among the transitional 
and quiescent galaxies, they still contribute at most 50% of 
Lir. The contribution of AGN Lir from these galaxies to the 
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Fig. 23. — Fraction of the dust luminosity due to ambient ISM as a function 
of rest-frame galaxy color Total dust luminosity includes includes energy 
absorbed in birth clouds (< 10 Myr) plus the ISM. The ratio is constrained 
by dust prescription of Chariot & Fall (2000). Even blue-sequence galaxies 
can have a high fraction of their dust luminosity absorbed away from birth 
clouds. 

global SFR density is beyond the scope of this work, but is 
most likely not very high. 

7. DISCUSSION 

7.1. Dust heating in actively star-forming galaxies 

Analysis presented in ^ indicates that the IR luminosity 
extrapolated from mid-IR flux is better correlated with the op- 
tical light of intermediate age populations than with the UV 
light of young stars. This result can be interpreted as the larger 
contribution of intermediate-age stars than of young stars in 
the mid-IR dust heating. Such interpretation is at odds with 
recent studies on nearby galaxies that find very good corre- 
lation between nebular emission line Paa, that comes from 
massive yo ung stars (< 10 Myr old), and the rest-frame 24^(m 
luminosity (ICalzetti et aP 120071: lAlonso-Herrero et al.ll2006c 
iRieke et al.ll2009ir Note, however, that we explore a differ- 
ent part of the mid-IR wavelength range (10-18 jim) which is 
more affected by PAH emission, and could therefore be more 
strongly correlated with the cold, diffuse dust from older stel- 
lar populations. What fraction of total Lir (irrespective of IR 
wavelength range) is expected to come from young stars in 
this particular sample? The dust model that we use in our 
SED fitting (Chai-lot & Fall 2000) allows us to estimate the 
relative contribution of the stellar energy absorbed by the stel- 
lar birth clouds (emitted by young stars) and the ambient ISM 
(emitted by intermediate age and old stars). This is achieved 
simply by considering the UV/optical luminosity that is ab- 
sorbed in these two components. In Figure |23] we plot the 
fraction of dust luminosity contributed by the ISM, i.e., away 
from the sites of current star formation. Not surprisingly, this 
fraction correlates well with NUV-/? color, which to first or- 
der gives the ratio of the recent to past star formation. For 
blue, star-forming galaxies ("(NUV-/?) < 3.5) the fraction 
of dust heating, and therefore the Lir due to ISM can be as 
high as 60% and is typically 40%. While still not dominant, 
the ratio of dust heating due to ISM can thus be quite signif- 
icant. Note that the ratio presented here does not constitute a 
measurement, but is set by the dust model we use here. How- 
ever, this dust model is physically motivated and can therefore 
serve as a guide. In reality, the contribution of the ISM may be 
somewhat different. Generally, there are many uncertainties 



with respect to the evolution of the intermediate age popula- 
tion and their dust production to leave room for their greater 
contribution to dust luminosity, perhaps even at mid-IR wave- 
lengths. 

7.2. Dust heating in quiescent galaxies 

The presence of high IR luminosities in galaxies that appear 
quiescent (not actively star-forming) based on UV/optical 
SEDs seems puzzling. Since for a given current SFR we 
have such a wide range of IR luminosities (including a large 
number of optically luminous, yet 24 /im-undetected galax- 
ies [Fig. |5l), the disconnect between the star formation and 
the IR properties appears strong. We have shown that AGNs, 
while contributing, do not dominate in the mid-IR, and that 
the absorbed stellar emission (from intermediate age and older 
stars) can on average reproduce high IR luminosities. Yet, this 
requires attenuations that are significantly higher than what 
we see in nearby galaxies with similarly low levels of spe- 
cific star formation. Locally, such galaxies are morpholog- 
ically early-type galaxies, with low to moderate amounts of 
dust and with infrared lumino sities that do not exceed 10 "'L© 
(iGoudfrooij & de Jongll995h . Cursory examination of optical 
HST ACS images of galaxies in our sample with log Lir > 10 
and "(NUV -R) > 5 (corresponding to UV/optical colors of 
nearby ellipticals, iDonas et aP 1200 7) indicates that 2/3 of 
them indeed look like early-type galaxies, and the rest are ei- 
ther edge-on disks, or show some structure. These results in- 
dicate that a fraction of early-type galaxies at higher redshifts 
have significantly higher dust contents, leading to higher IR 
luminosities. The presence of large amounts of dust even in 
nearby ellipticals is an open question (Temi et al. 2007), and 
is outside of the scope of this work. 

Another explanation for apparent high IR luminosities is 
that because we use IR templates based on star-forming galax- 
ies to estimate the total Lir of more quiescent galaxies, that 
this leads to significant overestimates. This explanation is not 
intuitive since one expects quiescent galaxies to have colder 
dust and therefore the mid-IR flux point used in conjunction 
with star-forming templates to underestimate the total Lir. 
But this explanation could be valid if quiescent galaxies con- 
tained significant contribution of stellar components that peak 
in the mid-IR (and are not incl uded in IR template s), such as 
the dust around the AGB stars (iBressan et al.lll998h . 

7.3. Monochromatic IR luminosity and SFR 

Recently there have been efforts to explore the use of 
monochromatic mid-IR luminosity as a tracer of star for- 
mation, either as a substitute for the total IR luminosity or 
as a measure that is intrinsically be tter correlated with SFR 
dCalzetti et alj2007l:lRieke et al.l200 9). Therefore, for our ob- 
servations at 24 /im, we construct a luminosity estimate at 12 
/im (L12), which corresponds to rest-frame of our observations 
at z = 1. To get L12 for g alaxies at other redshifts, we again 
rely on IR templates (Dal e & Heloull2002h . but now to obtain 
only a relatively small K-correction, instead of a full bolo- 
metric correction. Comparing L12 to FUV and B-band dust- 
corrected luminosities we find that the linear fits have scatters 
of 0.34 and 0.28 dex, respectively, i.e., they are some ^ 15% 
smaller than in relations with Lir. However, this comparison 
can be misleading since the range of Ln values is different 
(smaller) than of Lir. If instead we compare Pearson correla- 
tion coefficients, we find that they are basically the same for 
L12 and for Lir. While this does not mean that the total lu- 
minosity is intrinsically better correlated with the UV/optical 
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Fig. 24. — Mid-IR-excess (IRX) calculated according to IDaddi et alj 
(2007b) as a function of IR luminosity. IRX is defined as the ratio of SFR 
summed from Lir plus uncorrected FUV SFR to dust-coiTecte d FUV SFR 
(see text for details). This plot should be compared to Fig. 2 of'E)addi et al.' 
|2p07b), which shows many z ~ 2 IR-excess objects (IRX> 3.16, above the 
dashed line) with ULIRG luminosities, while there are few such objects in 
this sample. Nevertheless, we find some objects with moderate IR excess. 
The large majority of AGNs (88%) are not IR-excess objects, although AGNs 
are somewhat more frequent among IR-excess objects (10%) than among 
those that are not (4%). 

SFR than the 12 fim luminosity, it does at least argue that at 
our level of precision the correction to total luminosity neither 
does introduce significant additional uncertainty, nor does it 
offer any measurable benefits. 

7.4. IR excess and Compton-thick AGN 

Recently. iDaddi et alJ ( |2007 b, D07) have studied a popula- 
tion of z ~ 2 galaxies in GOODS fields that exhibits a mid-lR 
excess (IRX) around 8 /imand ascribed this excess to heating 
from Compton-thick AGNs. In ^we argued against the need 
for non-stellar sources of IR heating, yet one would like to ex- 
plore if there is a population of similar mid-lR-excess sources 
in our sample (at 10-18 /Ltm). D07 define mid-lR-excess ob- 
jects as those with the ratio of the combined IR and UV SFR 
(the latter not corrected for dust) to dust-corrected UV SFR 
exceeding 3.16 (10° ^): 

IRX(D07) = SFRiR+uv/SFRuv.corr > 3.16. (5) 

D07 obtain IR SFR by extrapolating 24 /im flux (which cor- 
responds to 8 /i m rest-frame flux ) to total IR luminosity 
and then using the iKennicutj (IT998h SFR conversion. They 
get dust-corrected UV SFR by applying a correction based on 
the fixed relation between the UV slope and UV attenuation. 
We will construct the IRX measure in exactly the same way, 
except that our Lir is extrapola ted from 10-18 nm rest-fr ame 
flux. Using models based on iBruzual & Charlo3 (l2003l) we 
find the following relation ("K-correction") between the ob- 
served B-z color at z = 2 and the rest- frame UV color: 

(B-z),=2= 1.8°(FUV-NUV). (6) 

Therefore their relation between the dust red dening and 
the ob served B and z magnitudes at z = 2 (iDaddi et al.l 
l2007ai Eqn. 8) corresponds to FUV attenuation of Apuv = 
4.5 (FUV— NUV), wh ich is similar in slope to that proposed 
by iMeurer et al.l (1199 9) for star-bursting galaxies. In Figure 
l24l we plot IRX against Lir for our galaxies. This figure 
should be qualitatively compared to D07 Fig. 2 except that 



the horizontal axis in D07 shows 8 /im rest-frame luminos- 
ity (L(8 /im)) while our figure shows total Lir. The Lir range 
of our sample of 9 < log Lir < 12 roughly translates into the 
range 8.5 < logL(8/im) < 11 range (Daddi et al. 2007a). The 
dashed line designates the IR excess criterion of D07. The 
vast majority of galaxies in our sample show no IR excess 
(IRX(D07)w 1). A small number of galaxies in our sample 
are found above the limit. Unlike in D07 where IR-excess ob- 
jects comprise nearly all high-luminosity objects, here there is 
almost no range in IR luminosity where excess objects domi- 
nate. Also, the maximum levels of IR excess in our sample are 
around 30, while in D07 they are up to 10 times higher This 
is probably related to the very different nature of the two sam- 
ples. Majority of the extreme IR-excess objects of D07 have 
logL(8/im) > 11, which corresponds to ULIRG luminosities. 
Such objects are all but absent in our sample. Indeed, a few 
objects that have Lir > 1O'^L0 in our sample have a mod- 
erate excess. Trying to determine if the IR-excess objects in 
our sample harbor Compton-thick AGN is beyond the scope 
of this paper Some recent work suggests alternative expla- 
nations for mid-IR excess at z ^ 2 involving PA H emission 
at 8 um dMurphv et alJl2008t iHuang et alJl2^ . The AGN 
that we identify in Section|6]are mostly not IR-excess sources 
(only 12% of optical or X-ray AGN have IRX> 3.16), al- 
though they do represent somewhat higher fraction among IR- 
excess sources (10%) than among those which are not (4%). 
This is consistent with conclusions of Section |6] that on av- 
erage AGNs contribute moderately to the mid-IR flux, unlike 
the possibly dominant AGN contribution among the IR lumi- 
nous mid-IR excess sources of D07. Also, note that Eqn.|5] by 
applying the Kennicutt conversion between IR luminosity and 
SFR, assumes that IR luminosity is dominated by young pop- 
ulations. This is certainly not true for more quiescent galaxies 
present in our sample. 

8. CONCLUSIONS 

In this paper we present (1) SFRs based on the Bayesian 
SED modeling of the UV and optical stellar continuum emis- 
sion, obtained by applying the Chariot & Fall (2000) dust at- 
tenuation model to a suite of Bruzual & Charloti (l2003i) stel- 
lar population synthesis models, and (2) the total infrared lu- 
minosities extrap olated from 24 /tm o bservations using the 
SED templates of iDale & Heloul (|2002|) . calibrated with local 
luminosity-color relations. We cover the redshift range up to 
z = 1.4, and study galaxies with spectroscopic redshifts from 
the AEGIS survey (Fig.[Tll. Our sample of 24 //m-detected R 
band-selected objects contains normal star forming galaxies 
and LIRGs, as well as quiescent galaxies, and is not biased 
against IR-luminous populations (Figs.|6]and[8]|. We compare 
IR luminosity with UV/optical SFRs averaged over various 
timescales, thus probing the present-day luminosities of stel- 
lar populations ranging in age from 0.1 to 3 Gyr. From this 
analysis we conclude the following: 

1 . When comparing UV/optical SFRs to IR luminosities, 
we confirm that one needs to treat actively star forming galax- 
ies separately from more quiescent ones. Th is caveat is wel l 
known in the study of nearby galaxies (e.g., Kennicuttll 19981) 
but is sometimes neglected at higher redshifts. Points 2-6 be- 
low pertain to actively SF galaxies which we select using a 
cut on rest-frame NUV-7? color (Figs.|5]and|6]. 

2. UV/optical SFRs averaged over relatively short 
timescales (10^ yr), and thus representing current SFRs, com- 
pare well (average difference 0.03 d ex) with Lir converted 
into SFR using the iKennicutI (119981) conversion. However, 
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the scatter between such SFRs and Lir is relatively high (0.42 
dex) (Fig.©. 

3. The scatter between UV/optical SFRs and Lir reduces 
as one considers SFRs averaged over longer periods of time, 
and is best for timescales between 1 and 3 Gyr, depending on 
the color (i.e., dominant population age) of a galaxy (Figs.fTol- 
[13] Equivalently, but more fundamentally, this means that the 
correlation is better between Lir and B-band dust-corrected 
luminosity (0.32 dex) than against FUV dust-corrected lumi- 
nosity (0.39 dex) (Fig. [I4]|. This argues for a significant role 
of intermediate-age stellar populations in mid-IR heating. 

4. Better correlation of Lir with optical luminosity than 
with FUV luminosity holds in redshift bins throughout 0.2 < 
z < 1.4, corresponding to 10-18 /im rest-frame wavelengths 
(Fig.IBJ. 

5. For our sample, which mostly consists of LIRGs and 
normal star forming galaxies, we find that galaxies are on 
average not optically thick, i.e., their IR luminosity (extrap- 
olated from the 24 /im flux) can be estimated from UV and 
optical photometry to within a factor of two (Fig.O . 

6. Many green valley galaxies are simply dust-obscured 
actively star-forming galaxies. However, there exist 24 /im- 
detected galaxies, some with LIRG-like luminosities, which 
have little current star formation (low specific SFR), i.e, they 
belong to green valley and even the red sequence because of 
their star formation history, not just dust reddening (Figs. [16] 
and[l7]l. 

7. On average, modeled amounts of dust absorption of 
stellar light are sufficient to produce the observed levels of 
IR luminosity, both for blue and for red-sequence galaxies 
(Fig.[T8]l. For red, quiescent galaxies this must include a large 
contribution of intermediate and old stellar populations and 
higher dust attenuations than in nearby early-type galaxies. 

8. Identified AGNs on average do not contribute signif- 
icantly to mid-IR luminosity at these redshifts. We see no 
evidence for a contribution by optical (type 2) AGNs to Lir 
and only up to 50% contribution by X-ray selected AGN, 
primarily at intermediate specific SFRs. Individual galaxies 
where AGN contribution to Lir is around 90% are not very 
common (Figs. [T9]j22] |. 

9. Extreme IR-excess sources similar to those identified 
at z ^ 2, and possibly related to Compton-thick AGNs, are 
very rare in our sample. Moderate IR excess can be attributed 
to either intermediate-age stellar populations or moderate IR 
heating from AGN (Fig.[24]i. 

10. Our findings (items 1-9) are qualitatively th e same 
if we compute IR luminosities using ICharv & Elba3 dlOOT]) 



IR templates instead of [Dale & Heloul ( l2002l) . However, the 
SED-derived SFRs becomes increasing ly degenerate with re- 
spect to Lir if the Lir is computed using iRiekeet all (l2009l) IR 
templ ates (which yi eld Lir up to an or der of magnitude hi gher 
flian lDale & Heloull2002i) . Thus, with RiekeetaD (IIQQI) de- 
rived Lir item 5 would no longer hold (Fig. [7]). 

1 1 . Using a fixed correlation between the UV color (spec- 
tral slope) and the FUV attenuation to obtain a dust correction, 
as opposed to a correlation that takes into account the effects 
of SF history on UV color, has the effect of producing dust- 
corrected SFR estimates that are on average better correlated 
with SFR averaged over 10*^ yr, than over 10^ yr, i.e., such 
procedure makes UV SFR a poorer indicator of the current 
SFR (Appendix B). 

Our work offers a new approach to study the relation be- 
tween star formation and infrared heating. The results would 
have implications for a number of studies which use mid-IR 
luminosity as a tracer of the current star formation. For ex- 
ample, it could affect the "time resolution" of cosmic SFR 
densities derived from 24 /im data. Our results are empiri- 
cal and are derived from typical data sets used at intermediate 
redshifts, but in future work we intend to extend this study to 
other redshift regimes and by employing other star formation 
indicators. 
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APPENDIX 

ROBUSTNESS OF STAR FORMATION RATES FROM SED FITTING 

Many of the results presented in this work depend on UV/optical star formation rates that we derive using the Bayesian SED 
fitting. In this section we evaluate the robustness of our SED fitting technique with respect to SFRs. We achieve this through 
simulations in which we try to recover a known SFR. In order to make the simulation as appropriate for our sample as possible, 
we proceed in the following manner Our SED fitting using real data tells us which model in our library best fits a given real 
galaxy. So in simulation, we simply substitute the observed fluxes with corresponding model fluxes, but using the observed flux 
errors and the redshift. Model fluxes are scaled to match the observed ones in / band. Then, the simulated SED fitting proceeds 
as it would for a real galaxy except that we exclude from the library the model whose fluxes we are trying to fit. In this way one 
gets an exact representation of the SED fitting for our sample but with the advantage that the true SFR (and any other parameter) 
that one tries to recover is actually known. 

First, we look at how well can the current SFR be recovered. By current, in this study we mean the star formation rate averaged 
over fsF =10^ yr, the shortest timescale corresponding to non-ionizing UV emission. In Figure[25]/e/f we show the comparison of 
the current SFR retrieved from SED fitting and the true SFR over the same timescale (fsF =10^ yr). The objects shown correspond 
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Fig. 25. — Comparison of SFRs recovered from UV/optical SED fitting vs. the true current SFRs. The true cun'ent SFR is known from the models. The SED 
SFR is averaged over 10^ yr (left) or over the population age (1-3 Gyr) (right). As exp ected, the cun'ent SFR is best recovered by a SFR averaged over a short 
timescale (left). If Lir followed the cuo'ent SFR, one would expect Fig. llOl and Fig. llll to respectively look like the left and the right panel here, flipped around 
the x = y axis. The fact that they show the opposite behavior argues that Lir does not correlate the best with the current SFR. 

to the sample of blue galaxies detected at 24 /im (i.e., the sample studied in most of ®. We recover the SFRs reasonably well. 
On average, the SED SFRs fall 0.04 dex below the true ones. The level of discrepancy depends on the SFR itself. For SFR 
between 1 and 10 Moyr"' there is on average no discrepancy, while for SFR between 10 and 100 MqJT'^ SFRs from SED fitting 
are 0. 13 dex lower than true ones. It should not be surprising that systematics at this level are present. The cause is most likely the 
limitations in the ability to obtain the full dynamic range of dust attenuations as discussed in S07. In any case, the main results in 
this work rely not on absolute rates of SF, but on their scatter vs. Lir. The scatter between SED SFRs averaged over 10^ yr and 
the true current SFRs is only 0.23 dex. This is very close to what is expected from our estimates of the individual errors of SED 
SFRs (Fig.EJi. 

Next we investigate if the derived SED SFR averaged over a timescale longer than 10^ yr can correlate better with the current 
SFR than the SED SFR averaged over 10** yr itself. This could possibly be the case if our SED SFRs over 10** yr are simply 
more noisy than SFRs averaged over longer timescales. One can imagine that this could result from the stochastic nature of SF 
histories in our models, where bursts have a timescale of roughly 10 yr. In Figure |25] ng/if we show how SED SFR averaged 
over the population age (1-3 Gyr) compares to the true current SFR. The scatter, 0.39 dex, is considerably worse than in the case 
where we averaged SED SFR over the short timescale. Indeed, by checking SED SFRs averaged over other timescales as well 
(1 and 2 Gyr), we find that the current rate of SF is indeed best recovered with the lO'^-yr timescale. Additionally, we also find 
(but do not show in a plot) that the true SFR averaged over the population age itself is very well recovered in SED fitting, with a 
negligible systematic offset and a scatter of only 0.20 dex. 

DUST CORRECTION OF UV/OPTICAL FLUX USING THE UV SLOPE 

To construct the dust-corrected SFRs and UV/optical luminosities we apply the attenuation model of ICharlot & Falll (l2000h 
directly to stellar population models and then compare the reddened models with the observed SEDs to derive SFRs and other 
parameters. Thus we use the full UV/optical SED to constrain the dust attenuation. This procedure is equ ivalent to an implicit use 
of the correlation between the UV slope (i.e., UV color) and the FUV attenuation (iCalzetti et al.llI994l) . but is not identical to it 
since the Chariot & Fall (2000) model (and implicitly our SED fitting) accounts for the effects of the galaxy SF history on the UV 
color (Buat 1992). However, in many instances it is more practical to perform the dust correction explicitly without considering 
the effects of SF history. In those cases one is using some^jcet/ correlation between the UV slope and the FUV attenuation. For 
our sample of blue-sequence (°(NUV-/?) < 3.5) galaxies we find that the best fixed-slope relation can be fit with 

Af,/y = 3.68''(FUV-NUV) + 0.29, °(FUV-NUV) < 1, (Bl) 

where ^{FUV - NUV) is the rest-fr ame UV color This relation i s somewhat steeper t han the equivalent relation for local (z ~ 0.1) 
SDSS galaxies (Salim et al.l l2007h . but still not as steep as the i Meurer et al.l d 19991) relation for local starbursts. While we use 
the exact same stellar population and dust models as in S07 and the fit is constructed in the same way (linear fit through running 
medians), the current color cut is somewhat bluer, and more importantly, the two samples are different, with local galaxies being 
more quiescent on average. 

Here we would like to draw attention to a systematic effect that, to our knowledge, has not been discussed elsewhere. Namely, 
while one normally expects the unattenuated FUV flux to best correlate with the SF on timescales of 10^ yr, using the fixed-slope 
relation to correct the UV flux (i.e., using the same relation between the UV spectral slope (or color) and the attenuation for 
all star forming galaxies, irrespective of their SF history) will effectively produce a measure of SF that is on average somewhat 
better correlated with the SF averaged over 10'' yr than, as expected, over 10^ yr For our sample of blue galaxies, we find the 
scatter around the linear least-square fit between SFRs (from SED fitting) averaged over 10^ yr and the fixed-slope dust-corrected 
FUV luminosity to be 0.15 dex, compared to 0.17 dex for SFR (from SED fitting) averaged over 10'^ yr On the other hand, as 
expected, the correlation of FUV luminosity corrected with full SED fitting dust-corrected FUV is the best (0.1 1 dex) for SFRs 
averaged over 10** yr, and significantly worse (0.19 dex) for 10^ yr The likely explanation for this counter-intuitive effect is that 
using a fixed slope between the attenuation and the UV color has the effect of overestimating the attenuation for galaxies that 
in reality lie below that fixed slope, and underestimating it for those that lie above it. Since galaxies below the fixed slope are 
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more likely to be galaxies with declining SF, while those above it tend to be more bursty (e.g. jKong et al.ll2004l) . the SFRs of the 
former get boosted (and thus become closer to an average over a longer period), while SFRs of the latter are suppressed, again 
mimicking the average that includes the pre-burst period. Because using the fixed slope effectively lengthens the SF timescale 
for galaxies with rising or falling SF histories, it has a consequence (given what we have shown in ® that it will correlate better 
with the IR luminosity than the FUV luminosity that was corrected using the full SED modeling. In other words, using a fixed 
slope to correct FUV dilutes our ability to constrain SF over various timescales and study its relation to the IR. 
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^ Galaxies matched to DEEP2 spectroscopic data set with (a) secure spectra and (b) lying in the intersection 
of CFHTLS and GALEX coverage (dark gray area in Figure [T). Note that a galaxy is kept in the sample 
regardless of the presence of a detection in a given UV or optical/near-IR band. 



